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Preface 


It  has  been  my  purpose  In  this  thesis  to  study  the  possibility  of 
finding  a  method  for  extending  spectrum  Identification  beyond  the  limi¬ 
tations  of  photopeak  analysis. 

To  carry  out  this  theoretical  study  the  Monte  Carlo  Transport  Code, 
Morse,  was  used  to  generate  output  spectra.  The  code  is  extremely  flexi¬ 
ble  and  has  been  adapted  for  Wrlght-Patterson's  CDC-6600  computer. 

I  cannot  express  my  gratitude  enough  to  the  personnel  of  RSIC  at 
Oak  Ridge  National  Laboratory.  Their  support  in  sending  me  Morse  and 
important  information  was  beyond  the  call  of  duty.  My  personal  thanks  to 
Dr.  Jabo  Sae  Tang  and  Dr,  M.  B.  Emmett  for  their  many  hours  of  getting 
me  through  the  troubled  days  when  I  could  not  get  Morse  running. 

A  little  closer  to  home,  my  appreciation  to  Dr.  George  John  and  Dr. 
Don  Shankland  for  their  ready  answers  to  my  detector  and  correlation 
theory  questions. 

Last,  I  wish  to  thank  Major  McKee  for  proposing  this  topic  and 
serving  as  my  Thesis  advisor.  And  I  must  thank  him  for  keeping  some  of 
my  grander  expectations  in  line  with  reality. 
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Abstract 


^‘Thls  report  shows  that  gamma  spectra  Identification  by  total 
flux  correlation  can  be  used  to  extend  identification  range  over 
photo  peak  methods.  Identification  was  based  on  two  decision  rules 
both  employing  cross-correlation  coefficients.  The  labgest  coeffic¬ 
ient  (first  decision  rule)  matched  the  unknown  spectra  with  the  cor¬ 
rect  source  thirty-seven  out  of  thirty-eight  trials.  The  proposed 
likelihood  function  (second  decision  rule)  had  a  success  rate  of 
thirty-five  out  of  thirty-eight  trials.  These  results  were  based  on 
spectra  generated  by  the  transport  code,  Morse. 

(• 
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I .  INTRODUCTION 


Background 

The  analysis  of  characteristic  gamma-ray  spectra,  from  spontaneous 
nuclear  decay,  for  identifying  the  atomic  source  or  sources  of  radia¬ 
tion  usually  employs  "matching"  the  unknovm  source's  spectral  features 
with  the  spectral  features  of  known  sources.  Figure  1  illustrates  this 
concept. 

Figure  1,  depicts  the  matching  of  the  unknown  source  in  a  rather 
simplified  manner.  The  match  shown  in  part  c  shows  the  spectra  being 
compared  by  their  energy  differences  of  the  peaks  and  relative  intensi¬ 
ties  of  the  peaks.  The  concept  of  the  library  and  measured  spectra  will 
be  expanded  upon  in  the  next  chapter. 

%  9  As  Illustrated  in  Figure  1,  the  spectral  features  that  are  most 

commonly  used  for  matching  are  the  peak  lines  in  the  spectra.  The  peaks 
are  the  easiest  means  for  several  reasons:  (1)  The  energy  of  the  peaks, 
or  differences  between  peaks,  can  be  determined;  and  (2)  the  relative 
intensities  of  the  peaks  can  be  determined. 

The  importance  of  knowing  the  relative  intensity  may  not  be  en¬ 
tirely  clear  at  this  point.  Basically,  what  if  two  sources  containing 
different  amounts  of  the  same  materials  were  measured?  A  proper  match 
would  then  depend  on  the  relative  peak  intensities,  as  well  as  line 
energies . 

As  a  review,  the  interactions  of  photons  with  matter  results  in 
either  partial  or  total  loss  of  energy.  The  three  major  processes  are 
(1)  the  photoelectric  effect,  (2)  compton  scattering  by  electrons  in  the 
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atoms  of  the  material,  and  (3)  pair  production.  Compton  scattering  is 
the  only  process  which  does  not  absorb  the  entire  energy  of  the  photon. 
These  processes  occur  in  both  the  transport  medium  and  detector  material 
The  peaks  in  a  spectrum,  like  those  in  Flgu’*e  1,  result  from  those 
gamma-rays  which  deposit  all  their  energy  only  within  the  detector's 
active  volume.  The  peaks  produced  by  these  uncollided  gamna-rays  are 
known  as  photopeaks.  Typically  these  peaks  dominate  the  spectral  fea¬ 
tures  and  are  used  for  source  identification.  The  uncollided  and  col¬ 
lided  photons  are  illustrated  in  Figure  2. 
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Figure  2.  Collided  and  uncollided  photons. 

As  the  photons  of  a  source  are  transported  to  further  and  further 
distances,  the  probability  that  a  photon  will  interact  in  the  transport 
medium  increases,  causing  the  ratio  of  collided  to  uncollided  photons 
to  continually  Increase.  Eventually,  there  would  be  so  few  uncolllded 
photons  reaching  the  detector  and  depositing  all  their  energy,  that  the 
identification  of  the  peaks  in  the  recorded  spectrum  would  no  longer  be 
practical  for  identification  purposes. 
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Now,  if  the  photopeaks  in  the  spectrum  were  the  only  information 
available  in  a  gamma-ray  spectrum,  then  isolation  of  these  peaks  become 
more  and  more  difficult  as  detection  range  Increases.  Eventually  the 
photopeaks  would  no  longer  be  Identifiable,  since  the  intensity  of  the 
photons  would  decrease  with  distance.  This  suggests  that  the  spectrum's 
information  vanishes  whenever  the  peaks  vanish. 

This  degradation  of  peak  lines  can  be  seen  in  figures  3  and  4. 
Figures  3  and  4  show  output  spectra  for  source  photons  with  energies  of 
1.0  Mev  and  .5  Mev,  respectively.  The  transport  of  photon  lines  must  be 
represented  by  the  transport  of  energy  groups  and  not  of  lines.  The 
group  structure  is  due  to  the  use  of  multigroup  cross-sections  that  use 
discrete  energy  groups.  The  1.0  Mev  photon  is  in  the  energy  bin  having 
limits  of  1.0-.80  Mev  while  the  .5  Mev  source  has  bln  limits  of  .45-. 5 
Mev.  The  energy  bins,  as  given  by  the  group  widths  of  the  crns-section 
library,  are  contained  in  Appendix  A. 

If  the  photopeaks  were  used  as  the  only  information  in  determining 
a  match,  identification  would  be  limited  by  the  ability  to  Isolate  and 
identify  the  photopeaks  of  the  measured  spectrum.  This  thesis  takes  the 
position  that  counting  the  uncollided  as  well  as  the  collided  photons 
might  add  "structure"  to  the  spectrum  unique  for  different  sources. 

This  added  "structure"  may  then  be  used  for  extending  the  identification 
range . 

Objectives 

The  basic  objective  of  this  thesis  is  to  model  and  study  the  remote 


detection  of  radioactive  materials  by  measuring  their  naturally  emitted 
gamma-rays,  obtaining  spectra,  and  then  employ  spectral  analysis  to 
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Figure  3.  Fluence  spectrum  for  a  1.  MeV  source  in  normal  density 
air  (1.298  g/1)  at  radii  of  20  meters  (.2  mfp) ,  166 
meters  (1.5  mfp),  and  AOO  meters  (3.0  mfp). 
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obtain  the  range  limitations  at  which  the  spectra's  sources  can  be  de¬ 


termined  . 

This  thesis  studies  the  hypothesis  that  the  total  spectrum,  consis¬ 
ting  of  both  uncollided  and  collided  photons,  may  be  used  as  a  means  for 
source  identification  and,  therefore,  extend  the  identification  range. 

In  other  words,  do  scattered  photons  contain  usable  Information  and  if 
so,  how  can  this  Information  be  extracted? 


Scope 

The  theoretical  study  lays  the  foundation  for  further  work  by  pre¬ 
senting  an  identification  technique,  based  upon  proposed  decision  rules, 
which  can  be  employed  to  select  a  source  from  its  spectrum.  This  study 
will  investigate  and  study  the  ability  of  the  identification  method  to 
make  the  proper  source  selection  as  the  photons  are  detected  at  increas¬ 
ing  distances.  In  addition  the  question  on  whether  a  "new"  source,  not 
Included  in  the  library,  has  been  detected  will  be  investigated.  The 
results  will  also  be  related  to  detection  time. 

To  reach  the  objectives,  two  basic  schemes  were  analyzed.  The  two 
scenarios  were; 

(1)  If  the  measured  spectrum's  source  is  known  to 
be  one  of  the  library  sources,  how  does  increasing 
the  detection  range  affect  the  ability  of  the  de¬ 
cision  rules  to  make  the  correct  source  Identifi¬ 
cation?  This  scenario  will  consist  of  both  a  "good" 
and  "poor"  measurement  for  the  measured  spectrum. 

The  meaning  of  "good"  and  "poor"  will  be  explained 
in  Chapter  II. 

(2)  If  a  measured  source  is  not  in  the  library, 
does  the  decision  rules  indicate  that  a  non¬ 
library  source  is  present?  Again  a  "good"  and 
"poor"  measured  spectrum  will  be  used. 


Limitations  of  Problem 

This  thesis  utilized  a  Monte  Carlo  transport  code  for  the  genera¬ 
tion  of  all  output  spectra  (Ref  1) .  Gamma-rays  were  transported  in  in¬ 
finite  (no  ground  or  structural  Interaction),  homogeneous  air  which  con¬ 
sisted  of  79%  ^2  and  21%  O2,  by  volume  with  a  density  of  1.298  g/1.  The 
modeling  of  the  photon  transport.  Imposed  the  following  conditions  on 
the  spectra  outputs: 

(1)  Sources  were  isotropic  point  sources. 

(2)  A  spherical  detector,  completely  surrounding 
the  point  source  was  used  to  estimate  fluence 
(See  Figure  5) .  This  assumption  is  necessary  to 
reduce  the  number  of  transported  photons  to  a 
reasonable  number.  The  characteristics  of  the 
transported  spectra  should  remain  valid. 

(3)  A  time  Independent  detection  was  modeled, 
with  the  number  of  source  photons  increasing  or 
decreasing  to  model  relative  detection  time.  In 
other  words,  the  fluence  is  estimated  not  the 
flux. 

(4)  Background  was  not  considered. 

(5)  A  realistic  detector’s  response  function 
was  not  modeled. 

Each  of  these  conditions  Imposed  restrictions  on  the  results  and 
ultimately  on  the  achievements  of  the  objectives.  A  brief  remark  on 
each  of  these  conditions  follows. 

Chosing  an  isotropic  point  source  to  model  the  configuration  of  the 
radioactive  material  was  not  a  bad  assumption.  If  source  dimensions 
were  no  greater  than  1  meter,  even  the  shortest  transport  range,  i.e. 

20  meters,  is  much  much  greater  than  the  source  dimensions. 

Condition  2  results  from  the  use  of  a  spherical  detection  surface 
surrounding  the  point  source.  Figure  5  Illustrates  this  configuration. 


Figure  5. 


Geometrical  arrangement  used  for  estimating 
fluence  of  a  point  source  with  the  Monte 
Carlo  Code,  Morse. 


This  setup  was  chosen  due  to  It's  ability  to  reduce  the  variance  of  the 
fluence  estimations  (Ref  2) .  This  fluence  estimation  technique  Is  one 
of  several  standard  methods  (Ref  2)  . 

A  time  Independent  detection  means  that  the  code  Is  not  modeling 
continuous  photon  emissions  from  the  source.  Rather,  the  total  Inten¬ 
sity  of  the  source  spectra  Is  determined  by  the  number  of  photons  that 
is  input  to  the  code  for  transport.  The  Intensity  of  the  bins  in  the 
source  spectra  la  determined  by  specifying  the  relative  amounts  of  the 
photons  for  each  energy  group. 

The  limitations  Imposed  by  conditions  4  and  5  contain  the  greatest 
deviations  from  a  real  detection. 

Condition  4  is  rather  straight  forward.  No  background  was  folded 
or  artificially  Incorporated  into  the  generated  spectra.  This  lack  of  a 
background  allows  fc r  a  perfect  environment  to  be  modeled.  Of  course, 
this  condition  is  never  obtained  in  a  real  detection  environment.  This 
thesis  assumes  total  background  elimination  Is  possible  while  at  the 
same  time  not  destroying  the  spectral  structure  created  by  the  collided 
photon  contributions.  Remember,  it  is  the  scattered  photon  contribu¬ 
tions  to  the  spectra  that  Is  being  investigated  as  a  possible  means  of 
extending  the  identification  range  over  the  photopeak  range. 

Condition  5  also  Incorporates  the  Ideal  detector  Into  the  transport 
model.  A  response  function  of  1  for  any  photon  crossing  the  detector 
surface  was  used.  This  Ideal  response  indicates  that  the  efficiency  for 
all  photons  of  any  energy  is  the  same.  This  condition  is  not  true  for  a 
real  detection  case. 
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Layout  of  Thesis 

Chapter  II  Introduces  the  sources  used  for  the  spectral  analysis 
and  discusses  how  they  were  constructed  and  transported. 

Chapter  III  Introduces  the  concept  of  spectral  analysis  by  cross- 
correlation  and  a  likelihood  function.  The  concept  of  a  library  and 
measured  spectra  are  defined  and  Illustrated. 

Chapter  IV  presents  the  results  of  the  spectral  analysis  and  Illus¬ 
trates  a  method  for  determining  minimum  detection  time. 

Chapter  V  follows  with  conclusions  and  recommendations. 

The  appendices  of  this  thesis  are  located  at  the  back  of  this  paper 


and  they  are  ordered  In  the  same  sequence  In  which  they  are  mentioned. 


II .  SOURCE  TRANSPORT 


Source  Considerations 


The  sources  considered  In  this  study  were  based  on  the  following 
assumption: 

Different  compositions  of  radioactive  nuclear  materials  emit 
gamma-radiation,  that  when  detected  close  to  their  surface, 
give  spectra  of  unique  characteristics  as  do  single  isotopes. 

Only  three  isotopes  were  considered:  U-235,  U-238,  and  Pu-239. 

The  source  spectra  were  artificially  constructed  from  combinations  of 

these  isotopes  based  upon  several  considerations: 

(1)  The  untransported  source  spectra  should  reflect,  to 
some  degree,  the  relative  amounts  of  the  combined  isotopes. 

(2)  The  source  spectra  should  have  some  major  lines  in 
common  but  with  different  relative  intensities.  This  con¬ 
dition  will  help  determine  the  ability  of  the  analysis 
method  to  distinguish  two  different  amounts  of  the  same 
material. 

(3)  Some  scattering  in  the  initial  photon  emissions  from 
the  isotope  sources  were  factored  into  the  spectra.  This 
spread  was  so  the  source  spectra  would  better  model  spec¬ 
tra  at  the  surface  of  the  material.  It  should  be  stated 
that  a  rigorous  down-scatter  analysis  was  not  attempted 
because  the  design  of  the  sources  was  deemed  unnecessary 
for  this  theoretical  study. 

Even  though  the  source  spectra  were  constructed  based  upon  these 
considerations,  the  details  of  the  source  spectra  are  somewhat  arbi¬ 
trary,  since  the  analysis  technique  is  being  evaluated,  not  the  detec¬ 
tion  of  any  specific  sources. 

The  major  energies  of  the  isotopes  are  tabulated  in  Table  I. 
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Table  I.  Major  gamma-ray  energies  of 

radioactive  material  (Ref  3:4) 


Energy 

Intensity 

Isotope 

(keV) 

l/(g-sec) 

U-235 

185.72 

4.3E4 

U-238  * 

1001.10 

1.0E2 

766.4 

3.9E1 

Pu-239 

129.28 

1.4E5 

413.69 

3.4E4 

*  These  energies  actually  arise  from 
Pa-234m  daughter  of  U-238. 

Figures  6  through  8  indicate  the  source  compositions  used  in  this 
thesis.  The  figures  are  labeled  source  A,  source  B,  and  source  C,  res¬ 
pectively,  and  will  be  referred  to  as  such  throughout  this  thesis. 

Table  II  gives  the  percentage  basis  of  the  source  spectra. 

While  the  energies  for  the  sources  were  based  upon  Table  I,  the 
figures  show  group  structure  instead  of  individual  line  energies.  Re¬ 
call  that  the  bin  limits  may  be  found  in  Appendix  A. 

All  spectra  were  graphed,  throughout  the  paper,  with  ordinate  units 
of  counts/keV/source  photon  vs.  abscissa  units  of  keV.  The  ordinate 
values  are  AttR^  times  fluence  for  each  energy  group.  Spectra  plots  were 
done  in  this  manner,  so  that  the  spectra  at  different  ranges  could  all 
be  represented  on  the  same  scale. 

Spectra  Generation 

All  of  the  output  spectra  were  artificially  produced  by  the  trans¬ 


port  code.  As  such,  estimates  of  the  group  fluence  values  (means)  could 
be  generated  to  the  statistical  degree  desired  by  running  the  code  as 


F’luence  (  Counts/Kevy  photon 

10“  ‘  '  10*'  10“  10' 


Source  A 


Figure  6.  Fluence  spectrum  of  Source  A. 
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Source  B 


Figure 
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Table  II.  Spectra  Intensities  of  Sources  on 
a  percentage  basis. 


Group 

(KeV) 

Source 

Percent  Intensity 
A  Source  B 

Source 

1300 

.0 

.18E-4 

.0 

1000 

.06974 

.0 

.2564 

800 

.0349 

.717E-5 

.1282 

700 

.3489 

.0 

.0256 

600 

.0 

.0 

.0 

512 

.2324 

.0 

.0 

500 

.0 

.0 

.0 

450 

.0177 

.1043 

.0 

400 

.0465 

.1043 

.0 

300 

.0116 

.0 

.0 

200 

.1164 

.7914 

.5129 

150 

.0698 

.0 

.0513 

100 

.0 

.0 

.0256 

75 

.0232 

.0 

.0 

60 

.0233 

.0 

.0 

45 

.0116 

.0 

.0 

30 

.0 

.0 

.0 

20 

.0 

.0 

.0 
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many  times  as  needed.  Recall  that  the  more  samples  that  are  used,  the 
better  the  true  mean  estimate  can  be  found,  i.e.  the  variance  of  the 
mean  is  proportional  to  1/N,  where  N  is  the  number  of  estimated  values 
for  the  quantity  of  interest  (for  example  f luence) .  With  this  in  mind 
the  output  spectra  were  grouped  into  two  categories. 

The  two  categories  were  the  library  spectra  and  the  measured  spec¬ 
tra.  These  categories  have  been  previously  mentioned  and  will  now  be 
fully  defined. 

The  library  spectra  were  generated  from  the  sources  in  a  manner 
that  allowed  for  obtaining  an  estimation  of  the  true  fluence  means. 

That  is  to  say,  the  library  spectra  were  generated  by  averaging  the 
fluences  of  many  runs .  This  resulted  in  a  standard  deviation  of  the 
mean  being  obtained  with  each  group  fluence  (Remember  the  energy  bln 
structure) .  Of  course,  if  the  sources  were  available  for  a  real  detec¬ 
tion,  then  careful  measurement  of  these  sources  under  controlled  condi¬ 
tions  would  be  carried  out.  This  is  what  the  library  spectra  represents. 

In  contrast,  the  measured  spectra  were  generated  so  as  to  represent 
a  measurement  under  less  than  optimum  conditions.  Hence,  their  fluences 
would  be  expected  to  be  randomly  distributed  about  the  mean  fluences  of 
the  library  spectra;  an  exact  match  was  not  likely. 

The  reasons  for  generating  spectra  at  different  ranges  will  be  ex¬ 
plained  shortly. 

The  measured  spectra  were  further  broken  up  into  two  types:  the 
"good"  and  "poor"  measured  spectra. 

The  terms  "good"  and  "poor"  were  used  only  to  indicate  that  the 
"good"  spectra  were  generated  with  a  larger  number  of  source  photons 
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than  the  "poor"  spectra,  thus  the  "good"  spectra  should  match  the  li¬ 
brary  spectra  more  closely.  The  term  "poor"  does  not  indicate  that  the 
generation  of  the  "poor"  spectra  resulted  from  computational  error. 

The  method  for  generating  these  measured  spectra  are  given  below 
under  the  appropriate  subheading. 

Library  Spectra.  The  library  spectra  consist  of  the  three  sources 
transported  to  8  different  distances:  20,  41,  83,  124.5,  166,  245,  333, 
and  390  meters.  Each  output  spectrum  was  determined  by  making  3  runs 
with  ten  batches/run  of  1000  source  photons/batch.  The  group  fluence 
values  for  the  3  runs  were  averaged  and  the  fractional  standard  devia¬ 
tions  were  determined  by  batch  statistics  as  given  by  the  following: 

Fluence  Average  (Multiple  runs) 

N  N 

^  E'l 

i=i  i=i 

N  =  number  of  runs 

=  number  of  batches  for  run  i 
=  average  value  of  interest  for  run  i 

F  =  averaged  value  of  Interest  (fluence) 


Fractional  Standard  Deviation  (Multiple  runs) 


f .s.d. 


Standard  Deviation  of  the  Mean  (Any  single  run) 
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where , 

N  =  number  of  batches  (note  the  change  from  multiple  runs) , 

n  =  total  number  of  Independent  histories  (function  of  the 
number  of  source  photons) , 

n^  =  number  of  Independent  histories  in  i-th  batch, 

=  accumulated  estimate  In  i-th  batch. 

N 

Note  that:  , 

"-E"! 

1=1 

N 


where  F 

i»  j 

batch. 


is  the  estimate  from  the  j-th  history  in  the  i-th 


N 


j=l 


*  F, 


(for  one  run) 


where  F  is  the  mean,  averaged  over  n  histories. 
The  fractional  standard  deviation  is 


f .s.d. 


These  statistics  are  computed  by  an  analysis  subroutine  within  the 
transport  code  (Ref  4) . 
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So,  if  the  desired  f.s.d.  is  not  low  enough  (f.s.d.  might  be  higher 

for  some  bins  than  others)  another  run  may  be  done  in  order  to  reduce 

the  f.s.d.  and  therefore,  reduce  a=  giving  a  better  estimate  for  the 

r 

fluence  means. 

Why  were  the  library  spectra  generated  for  different  distances? 
These  range  generations  were  done  so  that  the  concept  of  using  total 
spectra  analysis,  i.e.  uncollided  and  collided  spectral  features,  could 
be  carried  out.  Figure  9  best  illustrates  this  concept. 

By  "transporting"  the  library  spectra  at  the  same  range  as  a  mea¬ 
sured  spectra  a  match  can  be  attempted  using  the  total  spectrum  with  the 
hope  for  extending  the  identification  range. 

The  purpose  of  the  library  spectra  was  to  act  as  the  spectra,  cor¬ 
responding  to  known  sources,  by  which  a  match  with  a  measured  spectrum 
would  be  attempted.  A  perfect  match  would  indicate  that  the  measured 
spectrum  was  from  the  same  source  that  generated  the  library  spectrum. 
But  this  perfect  match  is  not  likely  and  only  a  best  match  would  be  ex¬ 
pected  . 

This  type  of  procedure  allows  for  the  investigation  of  the  hypothe¬ 
sis  that  the  scattered  photons  contain  useful  information. 

Measured  Spectra.  The  measured  spectra  consist  of  the  sources 
transported  to  the  same  distances  as  the  library,  but  with  fewer  source 
photons.  The  purpose  of  the  measured  spectra  was  to  act  as  the  unknown 
sources  for  a  match  with  the  library  spectra.  The  measured  spectra’s 
sources  will  be  referred  to  as  unknown  sources.  Of  course,  the  unknown 
sources  are  known  to  us,  so  that  conclusions  on  spectral  analysis  can 


be  evaluated. 


Also  included  was  Source  X,  which  was  a  source  not  used  in  the  li¬ 
brary.  Source  X  was  used  so  scenario  (2)  (Chapter  I)  could  be  answered. 
Source  X  is  shown  in  Figure  10 . 

Two  sets  of  measured  spectra  were  generated;  one  referred  to  as  the 
"good"  spectra  and  the  other  referred  to  as  the  "poor"  spectra.  Only 
one  run  with  one  batch  was  used  in  their  generation.  Therefore,  no  es¬ 
timate  of  the  mean  for  each  group  fluence  value  was  possible.  Vlhat 
these  measured  spectra  gave  were  fluence  values  randomly  distributed 
about  the  means  of  the  library  values. 

"Good"  Spectra.  These  measured  spectra  were  generated  by 
using  1000  photons  as  the  makeup  of  the  source's  Intensity. 

"Poor"  Spectra.  These  measured  spectra  were  generated  by 
using  100  photons  as  the  makeup  of  the  source's  intensity. 

It  should  not  be  surprising  to  expect  that  the  "good"  spectra  will 
be  distributed  with  less  variance  about  the  means  of  the  library  (of  the 
same  source)  than  would  the  "poor"  spectra.  Table  III  clearly  shows 
that  this  statement  is  correct.  Table  III  gives  the  number  of  standard 
deviations  away  the  measured  fluences  of  source  A  are  from  the  library 
fluences  of  source  A  at  the  range  of  20  meters.  Clearly,  Table  III 
(For  the  complete  listing  see  Appendix  C)  shows  that  the  standard  de¬ 
viations  of  the  "poor"  spectra  are  larger  than  the  "good"  spectra,  on 
the  average.  As  range  Increases  the  differences  between  the  "good"  and 
"poor"  spectra  become  even  more  apparent. 

The  purpose  for  generating  two  different  measured  spectra  for  each 
source  was  to  find  out  what  the  analysis  procedure  reveals  when  the 
search  spectrum  consist  of  either  a  good  measurement  (good  statistics) 
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Source 


Figure  10.  Fluence  spectrum  of  I 


2A 


k 


Table  III,  Distribution  of  "Good"  and  "Poor"  fluences 
about  the  means  of  their  library  fluences. 
Source  is  Source  A  at  the  range  of  20 
meters.  Values  are  in  Counts/Photon. 
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Group  (KeV) 

Library  A 
(s.  dev)* 

MA  (Good) 
(d.f .m.)** 

MA(Poor) 
(d.f .m.)** 

1300 

.0 

.0 

.0 

(.0) 

n/a 

n/a 

1000 

5.94E-2 

6.195E-2 

4.66E-2 

(7.47E-3) 

(-.34) 

(1.713) 

800 

3.043E-2 

3.548E-2 

1.541E-2 

(6.24E-3) 

(-.81) 

(2.406) 

700 

2.947E-1 

2.894E-1 

2.528E-1 

(1.77E-2) 

(.301) 

(2.369) 

600 

1.161E-2 

9.587E-3 

0.0 

(5.98-3) 

(.339) 

n/a 

512 

1.811E-1 

1.808E-1 

1.539E-1 

(8.10E-2) 

(3.3E-3) 

(.335) 

510 

l,866E-2 

1.551E-2 

0.0 

(6.23E-3) 

(.5059) 

n/a 

450 

2.472E-2 

1.635E-2 

5.716E-3 

(7.87E-3) 

(1.064) 

(2.415) 

400 

7.068E-2 

6.534E-2 

1.087E-1 

(1.37E-2) 

(.0389) 

(-2.77) 

300 

6.154E-2 

5.293E-2 

7.729E-3 

(1.512E-2) 

(.5696) 

(3.56) 

200 

1.183E-1 

1.404E-1 

1.491E-1 

(1.764E-2) 

(-1.25) 

(-1.74) 

150 

9.528E-2 

7.627E-2 

1.094E-1 

(1.63E-2) 

(1.162) 

(-.859) 

100 

1.492E-2 

3.563E-2 

2.235E-2 

(8.72E-3) 

(-2.37) 

(.8510) 

75 

2.083E-2 

1.439E-2 

1.309E-2 

(5.46E-3) 

(1.180) 

(1.418) 

60 

2.531E-2 

1.687E-2 

3.929E-2 

(7.33E-3) 

(1.151) 

(-1.90) 

45 

1.524E-2 

6.451E-3 

3.507E-3 

(6.37E-3) 

(1.378) 

(1.839) 

30 

3.897E-4 

0.0 

0.0 

(9.048E-4) 

n/a 

n/a 

20 

0.0 

0.0 

0.0 

n/a 

n/a 

n/a 

10 

0.0 

0.0 

0.0 

n/a 

n/a 

n/a 

*  (standard  deviation  of  library  (counts/photon)) 

**  (deviations  from  mean  of  library) 
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or  poor  measurement  (bad  statistics) . 

The  library  spectra  are  shown  in  Figures  11  through  13  at  various 
distances. 

The  "good"  and  "poor"  spectra  can  be  found  in  Figures  14  through 


21. 
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Figure  14.  "Good"  spectra  of  Source  A  at  ranges  of 
and  330  meters. 
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Figure  18.  "Good"  spectra  of  Source  C  at  ranges  of  20 
and  330  meters. 
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Figure  21.  "Poor"  spectra  of  Source  X  at  ranges  of  20,  166, 
and  330  meters. 


III.  CROSS-CORRELATION 


Background 

Cross-correlation  may  be  defined  as  a  mathematical  technique  for 
determining  the  degree  of  similarities  which  exist  between  two  spectra. 
With  the  introduction  of  digital  computers,  cross-correlation  became  a 
common  means  of  comparing  a  measured  spectrum  with  a  set  of  library 
spectra  in  search  of  sought-for  spectral  features.  The  main  application 
of  cross-correlation  techniques  has  long  been  to  detect  signals.  The 
first  discipline  to  widely  use  correlation,  and  probably  the  field  that 
still  uses  it  most  today,  is  the  communications  field.  The  application 
of  recognizing  periodic  communication  type  signals  in  noise  was  discus¬ 
sed  as  early  as  1950  (Ref  5) .  Several  books  that  deal  with  correlation 
techniques  have  been  written  with  application  primarily  to  communica¬ 
tions  (Refs  6; 7). 

Several  articles  were  found  that  specifically  applied  a  cross- 
correlation  technique  in  analysis  of  nuclear  spectra  (Refs  8;9;10;11). 
The  articles  presented  the  application  of  correlation  technique  in  iso¬ 
lating  and  enhancing  structure  in  pulse-height  spectrum  data.  The  spec¬ 
tra  peaks  were  the  main  structure  of  concern.  This  type  of  correlation 
method  is  a  peak  extraction  method.  A  peak  extraction  method  is  any 
technique  which  can  determine  if  a  gausslan  pulse  exists  within  the  ran¬ 
dom  noise  and  background  of  a  spectrum. 

An  article  in  1973  by  Horlick  (Ref  12)  presented  typical  results 
that  are  obtained  when  entire  spectra  are  cross-correlated  and  how  the 
I  resulting  cross-correlation  patterns  can  be  interpreted.  The  spectra 


used  were  line  emission  spectra  for  Co,  Nl,  and  Fe.  This  article  gave 
the  clearest  meaning  to  generating  library  spectra  of  knoim  sources, 
making  a  measurement  on  an  unknown  source,  and  then  matching  the  unknown 
spectrum  with  the  proper  library  spectrum.  The  match  then  Identifies 
the  source  of  the  unknown  spectrum. 


Decision  Rules 

In  order  to  analyze  the  unknown  spectrum  and  make  a  choice  as  to 
Its  source,  some  form  of  a  decision  must  be  made.  Two  decisions  rules 
were  proposed  and  studied.  Both  of  these  decision  rules  employ  cross- 
correlation  coefficients. 

Cross-Correlation  Coefficient.  The  cross-correlation  coefficient 
Is  defined  here  as 


P 


G 


* 


(1) 


where, 


.  =  cross-correlation  coefficient  of  spectrum  1  correlated 

with  spectrum  j. 

G  =  number  of  energy  groups  being  used  to  define  the  correlated 
spectra. 

Cj^(l)  =  spectrum  I's  fluence  estimate  of  energy  group  k. 

Ck(j)  =  spectrum  J's  fluence  estimate  of  energy  group  k. 


P 


-I  A 


may  be  normalized  by  the  denominator 


and  for  a  spectrum  correlated  with  itself 


The  actual  definition  of  (1)  is  quite  flexible  and  may  be  seen 
written  in  different  forms  in  other  literature  (Ref  6) .  In  determining 
the  form  of  the  correlation  equation  a  simple  expression  was  desired 
that  could  be  applied  to  the  group  nature  of  the  spectra  as  well  as  be 
easily  analyzed  in  a  statistical  fashion.  Equation  (1)  meets  these  de¬ 
sired  conditions  for  the  coefficient  equation. 

The  value  of  p.  .  can  be  thought  of  as  a  measure  of  all  spectral 
J 

features  shared  by  the  two  spectra  on  a  percentage  basis.  The  closer 
j  is  to  1  the  more  features  the  two  correlated  spectra  have  in  common 
with  one  another. 

Largest  Coefficient.  In  most  applications  of  the  cross-correlation 
coefficient,  the  match  for  a  measured  spectrum  has  depended  upon  cross- 
correlating  the  measured  spectrum  with  the  library  spectra  and  choosing 
the  largest  correlation  coefficient  as  the  indication  for  the  correct 
match.  This  method,  choosing  the  largest  coefficient,  constitutes  the 
first  decision  rule . 

Before  the  second  decision  rule  is  presented  it  is  instructive  to 
introduce  the  meaning  of  the  library  matrix  and  the  measured  matrix. 

The  library  matrix  contains  the  coefficient  values  of  the  library 
spectra  correlated  with  one  another.  Since  three  sources  comprise  the 
entire  library,  a  matrix  with  dimensions  of  3X3  will  be  formed.  The 


library  matrix  for  the  3  sources  looks  like  the  following; 


LA,LA 

*^LA,LB 

^LA,LC' 

LB,  LA 

‘^LB,LB 

‘’lb,lc 

'lc.la 

^LC,LB 

^LC,LC 

i  I 


Pt *  T D  would  be  the  cross-correlation  coefficient  of  the  library 

LAi  Lo 

spectrum  A  with  the  library  spectrum  B. 

Each  row  vector  of  the  library  matrix  may  be  viewed  as  the  "compo¬ 
nents"  of  library  spectrum  A  (row  1),  library  spectrum  B  (row  2),  and 
library  spectrum  C  (row  3) . 

Due  to  the  nature  of  the  correlation  equation,  this  matrix  is  sym¬ 
metric,  i.e.  *  Plc^LA*  PlB,LC  =  ^LC,LB' 

diagonal  terms  will  have  a  value  of  1,  since  cross-correlating  a  spec¬ 
trum  with  itself  must  result  in  a  perfect  match  (See  page  for  norma¬ 
lization  procedure) . 

These  coefficients  each  have  a  standard  deviation  of  the  mean,  a-, 
due  to  the  nature  of  their  generation.  Note,  the  diagonal  a*s  are  equal 
to  zero  because  the  normalization  procedure  always  gives  a  value  of  one 
for  the  diagonal  coefficients. 

The  completely  described  library  matrix  looks  like  the  following: 


^LA,LA  -  '^LA,LA  ^LA,LB  -  '^LA,LB  '’lA,LC  -  °LA,LC 

^LB,LA  -  '^LB.LA  ^LB,LB  -  '^LB,LB  ‘^LB,LC  -  °LC,LC 

PrriTA'^'PTnTA  PTr’TD"*'*^ir’TD 


'LC,LA  -  '^LC,LA  ^LC,LB  -  ‘^LC,LB  ^LC.LC  -  '^LC,LC 


t 


As  with  the  coefficients,  the  standard  deviations  are  also  sjnnmet- 

ric,  i.e  “  °LB,LA’  '^LB.LC  “  °LC,LB’  °LA,LC  “  '^LC.LA' 

A  library  matrix  was  calculated  for  each  distance  the  three  sources 

were  transported. 

The  measured  matrix  contains  the  coefficient  values  of  the  measured 
spectrum  cross-correlated  with  the  library  spectra  at  corresponding 
transport  distances.  In  this  case,  a  matrix  with  dimensions  1X3  results, 
giving 

[^M"i",LA  ^M"i",LB  ‘^M’'i",Lcj 

where  is  read  the  cross-correlation  coefficient  of  measured 

spectrum  "i"  ("i"  indicates  the  correct  source  for  the  measured  spectrum) 
correlated  to  library  spectrum  A  and  so  forth. 

Likelihood  Function.  The  likelihood  function  is  a  probability  den¬ 
sity  that  will  be  used  for  the  second  decision  rule. 

The  likelihood  function  used  in  this  study  is  given  by 


P(M=i)  = 


_  1  /'^Mi,LA~^Li,LA  ^  /pMi,LB~PLi.LB\  ^  /pMi,LC~PLl.Lc\ 

^  \  '^Ml.LA  /  \  ‘^Mi,LB  /  \  '^Mi.LC  / 


'^Mi.LA  ‘^Ml.LB  °Mi,LC 


where. 


M  -  assumed  source  for  the  measured  spectrum 

P(M=i)  =  value  of  the  density  function  assuming  M's  source  is 
source  i. 


^Mi  Lj  “  standard  deviation  of  the  correlation  coefficient  of  the 
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assumed  source,  M,  about  the  library  spectrum  of  j  (See  Appen¬ 
dix  B  for  derivation) .  "i"  indicates  the  assumed  source  of  the 

unknown  spectrum. 

Clearly,  (1)  consist  of  three  gausslan  density  functions  (gdf)  mul¬ 
tiplied  together.  The  use  of  the  gdf  is  based  on  the  assumption  that 
the  normalized  coefficients  follow  a  gaussian  distribution.  This  assump 
tion  must  be  kept  in  mind  when  the  results  of  this  decision  rule  are 
presented . 

The  likelihood  function  does  several  things  that  the  measured  coef¬ 
ficients  (alone)  do  not.  The  function: 

(1)  allows  for  the  similarities  of  the  library  spectra  with  each 
other  to  be  mathematically  taken  into  account;  and 

(2)  takes  into  account  any  error  associated  with  the  library  spec¬ 
tra  and  the  assumed  source. 

The  use  of  the  likelihood  function  is  a  more  mathematically  correct 
method  for  determining  the  correct  match  through  cross-correlation  analy 
sis  (Ref  13). 
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IV.  RESULTS  AND  DISCUSSION 


To  demonstrate  the  feasibility  of  extending  spectrum  identifica¬ 
tion  by  correlating  spectra  consisting  of  uncollided  as  well  as  collided 
photons,  the  measured  spectra  (including  Source  X)  were  correlated  to 
each  of  the  library  spectra.  The  decision  rules  were  then  applied  to 
the  coefficients  to  see  if  the  correct  match  was  indicated. 

The  correct  source  of  the  measured  spectra  will  be  designated  with 
an  M,  for  measured,  followed  by  the  letter  indicating  its  source,  i.e. 

MA  means  the  measured  spectrum  of  source  A.  In  addition  a  numerical 
value  may  follow  M"i"  indicating  results  for  the  "good"  or  "poor"  mea¬ 
sured  spectra,  i.e.  MAIOOO  would  be  for  "good"  while  MAIOO  would  be  for 
"poor"  spectra.  This  nomenclature  Indicates  the  number  of  source  pho¬ 
tons  used  for  their  spectra  generation. 

The  results  of  the  decision  rules  are  presented  in  the  sequence  in 
which  the  scenarios  were  presented  in  Chapter  II. 

Coefficient  Matrices 

Tables  IV,  V,  and  VI  give  the  coefficient  values  for  the  library 
matrix,  the  "good"  measured  matrix,  and  the  "poor"  measured  matrix,  re¬ 
spectively.  All  the  coefficients  shown  are  normalized  coefficients. 

The  associated  standard  deviations  of  the  means  for  the  library 
coefficients  are  tabulated  in  Table  VII. 

Library  Matrix.  Let  us  first  look  at  the  library  matrix  (Table  VI) 
and  understand  some  of  the  trends  that  are  occuring.  The  most  noticable 
feature  about  the  coefficients  is  that  their  values  are  becoming  closer 
and  closer  together  as  range  Increases.  This  phenomenon  is  occurring 


Table  IV.  Coefficient  values  of  the  library 


matrix 

SOURCE 

1.0000 

.1478 

.4325 

.1478 

1.0000 

.8530 

.4325 

.8530 

1.0000 

RANGE  = 

20  meters 

1.0000 

.3814 

.4556 

.3814 

1.0000 

.8587 

.4556 

.8587 

1.0000 

RANGE  = 

41  meters 

1.0000 

.5031 

.5727 

.5031 

1.0000 

.8821 

.5727 

.8821 

1.0000 

RANGE  = 

83  meters 

1.0000 

.7079 

.7578 

.7079 

1.0000 

.9174 

.7578 

.9174 

1.0000 

RANGE  = 

124.5  meters 

LEGEND 

1.0000 

.8188 

.8591 

.8188 

1.0000 

.9469 

^LA,LA 

‘’la,  LB 

‘’la.lc 

.8591 

.9469 

1.0000 

RANGE  = 

166.  meters 

^LB.LA 

‘’lb,  LB 

‘’lb.lc 

1.0000 

.8735 

.9163 

^LC,LA 

‘’lb,lc 

‘’lc,lc 

.8735 

1.0000 

.9576 

.9163 

.9576 

1.0000 

RANGE  = 

246.  meters 

1.0000 

.9090 

.9547 

.9090 

1.0000 

,9626 

.9547 

.9626 

1,0000 

RANGE  = 

330.  meters 

1.0000 

.9138 

.9696 

.9138 

1.0000 

.9618 

.9696 

.9618 

1.0000 

RANGE  = 

390.  meters 

1.0000 

.9307 

.9810 

.9307 

1,0000 

.9559 

.9810 

,9559 

1.0000 
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Table  V.  Coefficients  of  "Good"  matrix 


B 

I 

B 


RANGE 


(meters) 

MAIOOO 

MBIOOO 

MClOOO 

20.0 

(.995, .422, .497) 

(.379, .999, .853) 

(.470,. 859,. 998) 

41.0 

(.995, .472, .539) 

(.496,. 999,. 883) 

(.592,. 860,. 997) 

83.0 

(.990,. 637,. 709) 

(.704,. 997,. 921) 

(.751,. 904,. 996) 

124.5 

(.986,. 759,. 804) 

(.818,. 997,. 942) 

(.853,. 961,. 997) 

166.0 

(.980,. 848,. 893) 

(.869, .992, .949) 

(.904,. 939,. 988) 

246.0 

(.985,. 884,. 939) 

(.901,. 984,. 918) 

(.943,. 950,. 982) 

330.0 

(.953,. 928,. 950) 

(.924, .990, .963) 

(.965,. 951, .977) 

390.0 

(.965,. 916,. 954) 

(.921,. 974, .939) 

(.964,. 948, .990) 

Table  VI.  Coefficients  of  "Poor"  matrix 

i 

RANGE 


(meters) 

MAIOO 

MBIOO 

MClOO 

20.0 

(.969, .500, .533) 

(.364, .997, .861) 

(.446,. 836,. 981) 

41.0 

(.971,. 549,. 600) 

(.490,. 941,. 866) 

(.530,. 927,. 985) 

83.0 

(.886, .753, .759) 

(.678, .967, .866) 

(.720,. 898,. 944) 

124.5 

(.907, .687, .735) 

(.797, .960, .897) 

(.710, .807, .831) 

166.0 

(.909,. 715,. 763) 

(.795,. 940,. 897) 

(.789,. 780,. 878) 

246.0 

(.859, .778, .842) 

(.722, .841, .789) 

(.919,. 913,. 944) 

330.0 

(.946,. 907,. 932) 

(.844, .961, .914) 

(.731, .696, .711) 

i: 
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Table  VII.  Standard  deviations  of  the  means  for  the 
coefficient  library  matrix 


I 

I 


LIBRARY 

SIGMA  RANGE 

=  20  meters 

A 

B 

C 

0.0 

8.525E-003 

7.873E-003 

8.525E-003 

0.0 

8.101E-003 

7.873E-003 

8.101E-003 

0.0 

LIBRARY 

SIGMA  RANGE 

=  41  meters 

A 

B 

C 

0.0 

1.013E-002 

1.002E.002 

1.0i3E-002 

0.0 

1.133E-002 

1.002E-002 

1.133E-002 

0.0 

LIBRARY 

SIGMA  RANGE 

=  83  meters 

A 

B 

C 

0.0 

1.262E-002 

1.336E-002 

1.262E-002 

0.0 

1.268E-002 

1.336E-002 

1.268E-002 

0.0 

LIBRARY 

SIGMA  RANGE 

=  124.5  meters 

LEGEND 

A 

B 

C 

- 

0.0 

1.940E-002 

1.800E-002 

'^LA.LA 

°LA,LB 

°LA,LB 

1.940E-002 

0.0 

1.957E-002 

1.800E-002 

1.957E-002 

0,0 

'^LB.LA 

'^LB,LB 

'’lb,lc 

LIBRARY 

SIGMA  RANGE 

=  160.  meters 

_ 

A 

B 

C 

LC^LA 

LC)LB 

LC  f  I#C 

0.0 

1.893E-002 

1.893E-002 

1.893E-002 

0.0 

1.865E-002 

1.893E-002 

1.865E-002 

0.0 

LIBRARY 

SIGMA  RANGE 

=  246.  meters 

A 

B 

C 

0.0 

2.315E-002 

2.414E-002 

2.315E-002 

0.0 

2.147E-002 

2.414E-002 

2.147E-002 

0.0 

• 

LIBRARY 

SIGMA  RANGE 

=  330.  meters 

A 

B 

C 

0.0 

2.769E-002 

2.863E-002 

2.769E-002 

0.0 

3.251E-002 

2.863E-002 

3.251E-002 

0.0 

LIBRARY 

SIGM(V  RANGE 

=  390.  meters 

A 

B 

C 

0.0 

3.371E-002 

3.166E-002 

3.371E-002 

0.0 

2.960E-002 

3.166E-002 

2.960E-002 

0.0 
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due  to  the  Increased  downscattering  of  the  source  photons  as  range  in¬ 
creases.  The  downscattering  photons  cause  the  group  bins  of  the  lower 
energies  to  fill  up  creating  a  higher  correlation  between  the  library 
spectra  due  to  the  spectral  differences  disappearing.  This  downscatter 
is  evident  in  any  of  the  spectra  figures. 

Measured  Matrix .  Tables  V  and  VI  give  the  measured  matrix  for  the 
"good"  and  "poor"  spectra,  respectively. 

Like  the  library  coefficients  the  measured  coefficients  are  becom¬ 
ing  closer  and  closer  as  range  Increases.  Of  course,  this  was  expected 
since  the  same  sources  are  being  used  and  the  spectral  features  of  these 
measured  spectra  are  disappearing  as  well. 

Library  vs.  Measured  Coefficients 

The  purpose  of  generating  coefficients  is  to  compare  the  library 
coefficients  with  the  measured  coefficients,  apply  the  decision  rules, 
and  study  the  indicated  match.  The  measured  coefficients  are  plotted 
along  with  the  correct  row  of  the  library  coefficients  in  Figures  22 
through  27  (Source  X  will  be  presented  later) . 

The  open  symbols,  connected  by  some  line  "type",  are  coefficients 
of  the  indicated  library  spectrum  correlated  with  itself  and  the  other 
library  spectra.  The  crossed-in  symbols  are  the  coefficients  of  the 
measured  spectrum  (unknown  source)  with  the  library  spectra.  If  the 
measured  spectrum's  source  is  the  same  assumed  library  source,  then  each 
paired  symbol,  one  open  and  one  crossed,  should  be  "close"  to  one  ano¬ 
ther  or  show  the  same  ordering  at  each  range.  These  figures  indicate 
that  these  trends  are  true. 

These  figures  show  several  things.  First,  the  "good"  spectra 


LEGEND 

♦P  ■=  p(MA,LA) 
a  p(MA,LB) 

^  -  p(MA,LC) 

_o_ii  p(la,la) 
-CL_-7  p(la,lb) 

-  p(LA,LC) 


1000 


?00.0 


300.0 


Figure  22.  Measured  matrix,  MA  (unknovm:  Source  A) 
library  coefficients  of  row  1. 
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I 
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Coefficients 


Range  (Meters) 

Figure  27.  Measured  matrix,  MClOO  (unknown:  Source  C)  vs 
library  coefficients  of  row  3. 


5A 


LEGEND 

©  =  p(MB,LA) 

(3  =  P(MB,LB) 

^  =  p(MB,LC) 

_0 _ =  p(LA,LA) 

p{LA,LB) 
O  —  p(LA,LC) 


-T - j - 

100.0  200.0  300.0 
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.  28.  Measured  matrix,  MB  (unknown:  Source  B)  vs.  library 
coefficients  of  row  1  (THE  SOURCE  A  ROW  IN  THE 
LIBRARY  MATRIX) 


400.0 
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coefficients,  at  each  range,  "closer"  than  the  "poor"  spectra  coeffi¬ 
cients.  This  difference  is  due  to  the  "poor"  fluences  being  randomly 
distributed  about  the  means  of  the  library  with  greater  deviation  than 
the  "good"  spectra.  This  distribution  of  the  "good"  vs  "poor"  has  al¬ 
ready  been  mentioned  in  a  previous  section.  Here  the  plots  illustrate 
the  difference  between  the  distributions. 

Another  feature  about  the  coefficients  is  that  as  range  Increases 
the  coefficients  begin  to  converge  on  one  another  making  it  difficult  to 
tell  the  coefficients  apart.  This  bunching  implies  that  identifying  the 
source  of  the  measured  spectra  will  also  become  more  difficult  as  range 
increases  and  spectral  features  disappear. 

Figure  28  shows  what  happens  when  the  measured  matrix  coefficients 
are  compared  to  the  wrong  row  of  the  library.  In  other  words  the  wrong 
source.  The  symbols  are  no  longer  paired  and  the  largest  measured  co¬ 
efficient  is  not  ...  Pvn,  is  the  coefficient  correspond- 

NIj  j  IjiIj  Ko  y  LiA  ^  L<15 

ing  to  the  assumed  Source  B  (wrong  in  this  case) . 

Results  of  the  Decision  Rules 

So  far  the  relationship  between  the  library  and  measured  coeffi¬ 
cients  have  been  shown.  The  decision  rules  are  now  applied  to  the  co¬ 
efficients  in  order  to  study  how  the  correct  match  is  affected  as  range 
increases.  The  maximum  studied  range  at  which  the  correct  match  was 
chosen  will  then  be  used  to  compare  the  photopeak  identification  method 
to  these  proposed  methods. 

Largest  Coefficient.  The  largest  coefficient  in  the  measured  ma¬ 
trix  will  determine  the  most  likely  source  for  the  measured  spectrum. 
This  decision  rule  has  been  used  in  previous  work  (Ref  12) .  By  looking 


Table  VIII  gives  the  value  of  the  likelihood  function  along  with 
the  maximum  the  function  could  possibly  be.  The  maximum  would  indicate 
that  all  coefficient  pairs  were  perfectly  matched,  i.e.  each  pair  is  ex¬ 
actly  equal  to  each  other,  and  therefore,  the  numerator  would  reduce  to 
one. 

By  noting  the  maximum  indicated  for  each  assumed  source,  it  can  be 
seen  that  the  values  of  the  likelihood  functions  are  weighted  different¬ 
ly  (due  to  the  process  used  in  determining  the  standard  deviations  used 
in  the  function) .  If  the  function  is  divided  by  its  maximum  value,  then 
a  percentage  will  be  found  that  indicates  the  degree  the  function  is  to 
a  perfect  match.  This  ratio  will  be  referred  to  as  the  normalized  like¬ 
lihood  function.  The  largest  normalized  likelihood  function  is  used  as 
the  second  decision  rule.  Table  IX  gives  the  results. 

Table  IX  indicates  the  wrong  matches  in  the  location  where  the  sym¬ 
bol  "+"  is  present. 

Largest  Coefficient  vs.  Likelihood  Function 

From  the  nrevious  two  sections,  the  identification  of  the  correct 
source,  by  either  decision  rule,  gives  comparable  results.  While  the 
rule  utilizing  the  likelihood  function  indicated  three  wrong  matches 
compared  to  one  wrong  match  for  the  coefficient  rule,  each  wrong  match 
was  indicated  only  for  the  "poor"  spectra.  In  general  terms,  the  two 
rules  must  be  considered  ..o  be  equal  in  their  success  for  choosing  the 
correct  source. 

Cross-Correlation  vs.  Photopeak  Identification 

In  order  to  make  the  comparison  between  the  cross-correlation 


Table  VIII 


(a)  Likelihood  function  values  for  measured  source  MAlOOO 


RANGE  ASSUMED 
meter  SOURCE 


/ASSUMEDX 
P* SOURCE  / 


/ASSUMED^ 
Pmax I  SOURCE  ) 


A 

1.4E+2 

3.6E+2 

B 

2.5E-369 

4.5E+3 

C 

1.5E-269 

2.7E+3 

A 

4.2E+2 

6.6E+2 

B 

2.5E-125 

1.4E+3 

C 

7.3E-78 

9.5E+2 

A 

1.2E+2 

3.2E+2 

B 

1.8E-36 

1.2E+3 

C 

8.4E-15 

4.2E+2 

A 

7,8E+1 

1 . 3E+2 

B 

2.2E-2 

1.4E+2 

C 

1.6E-2 

2.1E+2 

A 

9.5E+1 

1.2E+2 

B 

2.9E+1 

1.9E+2 

C 

6.9E+1 

1.7E+2 

A 

5.6E+1 

5.9E+1 

B 

2.8E+1 

1.6E+2 

C 

6.5E+0 

8.7E+0 

A 

4 . 2E+1 

4.6E+1 

B 

2.8E+1 

3.5E+1 

C 

2.9E+1 

3.2E+1 

A 

2.0E+1 

2.1E+1 

B 

2.7E+1 

3.5E+1 

C 

4.3E+1 

5.1E+1 
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Table  Vlll 


f  (c)  Likelihood  function  values  for  measured  source  MClOOO. 


I 

I 

I 


i  «i  * 

i 


RANGE 

ASSUMED 

/ASSUMED\ 

/ASSUMED  \ 

P 

meter 

SOURCE 

P\ SOURCE  / 

PmaxlsOURCE  / 

Pmax 

20 

A 

4.0E-71 

3.6E+2 

0. 

20 

B 

3.3E-13 

4.5E+3 

0. 

20 

C 

2.0E+3 

2.7E+3 

.752 

41 

A 

1.8E-48 

6.6E+2 

0. 

41 

B 

3.7E-4 

1.4E+3 

0. 

41 

C 

6.9E+2 

9.5E+2 

.729 

83 

A 

1.5E-8 

3.2E+2 

0. 

83 

B 

5.3E-K) 

1.2E+3 

.005 

83 

C 

4 . 2E+2 

4.2E+2 

.963 

124.5 

A 

9.2E-1 

1.3E+2 

0. 

124.5 

B 

9.0E+1 

1.4E+2 

.645 

124.5 

C 

2.0E+2 

2.1E+2 

.981 

166.0 

A 

2.7E+1 

1.2E+2 

.226 

166.0 

B 

l.lE+2 

1.9E+2 

.580 

166.0 

C 

1.6E+2 

1.7E+2 

.936 

246.0 

A 

4.4E+1 

5.9E+1 

.755 

246.0 

B 

l.lE+2 

1.6E+2 

.705 

246.0 

C 

8.5E+0 

8.7E+0 

.971 

330.0 

A 

4.1E+1 

4.6E+1 

.898 

330.0 

B 

2.9E+1 

3.5E+1 

.834 

330.0 

C 

3.4E+1 

3.2E+1 

.970 

390.0 

A 

2.0E+1 

2.1E+1 

.957 

390.0 

B 

3.0E+1 

3.5E+1 

.853 

390.0 

C 

5.0E+1 

5.1E+1 

.978 
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(d)  Likelihood  function  values  for  measured  source  MAIOO. 


RANGE 

ASSUMED 

/ASSUMED\ 

/ASSUMED \ 

P 

meter 

SOURCE 

P\ SOURCE  / 

PmaxVsOURCE  / 

Pmax 

20 

A 

1.8E+1 

l.lE+2 

.159 

20 

B 

5.8E-111 

8.8E+2 

0. 

20 

C 

4.8E-72 

5.1E+2 

0. 

41 

A 

1.4E+2 

2.1E+2 

.643 

41 

B 

1.5E-34 

2.7E+2 

0. 

41 

C 

8.0E-20 

1.8E+2 

0. 

83 

A 

l.OE+1 

l.OE+2 

.097 

83 

B 

3.2E-4 

2.2E+2 

0. 

83 

C 

3.0E-1 

8.0E+1 

0. 

124.5 

A 

8.3E+0 

4.1E+1 

.202 

124.5 

B 

6.9E-1 

2.6E+1 

.027 

124.5 

C 

4.4E-1 

3.9E+1 

.011 

166.0 

A 

6.1E+0 

3.9E+1 

.159 

166.0 

B 

9.1E-1 

3.6E+1 

.026 

166.0 

C 

l.lE+0 

3.2E+1 

.035 

246.0 

A 

4.3E+0 

1.9E+1 

.230 

246.0 

B 

4.8E-K) 

2.9E+1 

.165 

246.0 

C 

4.3E-K) 

1.6E+1 

.266 

330.0 

A 

1.3E+1 

1.5E+1 

.879 

330.0 

B 

5.6E-K) 

6.5E+0 

.865 

330.0 

C 

5.4E-K) 

5.9E+0 

.922 
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Table  VIII 


(e)  Likelihood  function  values  for  measured  source  MBIOO 


RANGE  ASSUMED 
meter  SOURCE 


/ASSUMED] 

pIsource  ) 


/ASSUMED\ 
Pmax\ SOURCE  / 


Table  IX.  Results  of  "matching"  using 
the  likelihood  function 


technique  and  the  photopeak  technique,  the  range  where  photopeak  iden¬ 
tification  can  no  longer  be  utilized  must  be  determined.  Some  subjec¬ 
tive  decision  had  to  be  made  based  on  the  spectra  generated.  By  looking 
at  the  "good"  and  "poor"  spectra,  the  range  of  166  meters  shows  that 
most  of  the  initial  source  structure  has  vanished.  A  way  to  interpret 
the  range  limit  for  photopeak  identification  is  to  picture  the  measured 
spectra  of  Source  A  and  Source  C  (Figures  6  &  8)  and  their  Source  Spectra 
(Source  B  is  excluded  because  of  its  lower  spectra  bins,  l.e.  easily 
identified  from  the  other  spectra  by  its  upper  bin  limit) .  At  this 
range  it  is  becoming  difficult  to  look  at  the  spectra  and  identify  the 
correct  source.  Therefore,  the  166  meter  range  will  be  considered  the 
maximum  range  limit  where  photopeak  identification  may  be  carried  out 
successfully. 

In  order  to  help  justify  this  range,  a  look  at  the  mean  free  paths 
for  a  couple  of  the  primary  energy  groups  .8-l,MeV  and  .15-.2MeV  trans¬ 
lated  into  1.5  and  2.8  mean  free  paths,  respectively.  This  distance 
does  not  seem  to  be  a  bad  choice,  since  the  probability  of  the  photons 
scattering  is  77.3%  and  93.9%,  respectively. 

The  decision  rule  results  indicate  a  rather  large  range  extension 
for  Identification.  This  extension  is  by  as  much  as  twice  the  photopeak 
range  limit,  i.e.  166  meters  compared  to  as  high  as  390  meters. 

Source  X:  The  Non-Library  Source 

Table  X  shows  the  measured  coefficients  of  measured  spectra  X  and 
XlOO  with  the  library  matrix.  These  coefficients  along  with  library  co 
efficients  are  plotted  in  Figures  29  through  34. 

Clearly  each  decision  rule  will  always  indicate  the  wrong  match. 


Table  X. 

Coefficients  of 

"Good 

matrix  of  source 

X 

RANGE 

(meters) 

MXIOOO 

20.0 

(.431, .624,. 649) 

(. 

41.0 

(.598, .669,. 701) 

(. 

83.0 

(.814,. 805,. 853) 

(. 

124.5 

(.899,. 848,. 890) 

(. 

166.0 

(.954, .895,. 934) 

(. 

246.0 

(.957,. 938, .967) 

(. 

330.0 

(.970, .935, .973) 

(. 

390.0 

(.965,. 916,. 954) 

"  and  "Poor" 


MXlOO 

238. . 164.. 557) 

536..  594.. 557) 

750..  860.. 903) 

888..  784.. 844) 

898. . 845.. 868) 

916..  839.. 894) 

976. . 923.. 977) 
N/A 
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I  I  I  I  I 


0.0  100.0  200.0  300.0  400 

Range  (Meters) 

Figure  29.  Measured  matrix,  MX  (unknown:  Source  X)  vs 
library  coefficients  of  row  1. 
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I  I  !  1  r 


0.0  100.0  200.0  300.0  400.0 

Ran|G;e  (Meters) 

Figure  30.  Measured  matrix,  MXlOO  (unknown:  Source  X)  vs 
library  coefficients  of  row  1. 
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£ 


CD® 


Coefficients 


0.0  100.0  200.0  300.0  400 

Ran^e  (Meters) 


Figure  32,  Measured  matrix,  MXlOO  (unknown:  Source  X)  vs 
library  coefficients  of  row  2. 
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Coefficients 


Figure  33.  Measured  matrix,  MX  (unknown:  Source  X)  vs 
library  coefficients  of  row  3. 


0.0  100.0  200.0  300.0  400.0 

Range  (Met ere) 

Figure  34.  Measured  matrix,  MXlOO  (unknown:  Source  X)  vs 
library  coefficients  of  row  3. 
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What  Is  needed  is  some  method  able  to  discriminate  Source  X  as  a  "new" 
source . 

Examining  the  coefficients  of  Table  X  shows  that  there  is  no  clear 
"dominant"  coefficient  at  the  low  range,  l.e.  less  than  124.5  meters. 

The  magnitudes  at  these  low  ranges  differ  from  the  large  values  found 
from  the  "good"  and  "poor"  coefficients  for  the  correct  source.  To  ar¬ 
gue  that  these  low  values  Indicate  a  "new"  source  would  be  hard  to  jus¬ 
tify.  Even  if  this  statement  could  be  justified,  this  low  range  is  with¬ 
in  the  photopeak  range  and  therefore,  adds  little  to  the  analysis  of  X 
spectra. 

Detection  Time 

The  minimum  count  time  based  on  a  spectrum  consisting  of  a  total  of 
1000  source  photons  will  be  determined.  The  count  time  determined  is  the 
minimum  since  the  same  results  are  desired  as  the  ones  with  the  "good" 
spectra.  The  1000  photons  were  chosen  due  to  the  simple  fact  that  1000 
source  photons  were  used  as  input  for  the  generation  of  the  measured 
spectra. 

The  minimum  count  time  will  be  based  on  a  g-sec  basis  so  that  know¬ 
ledge  of  the  source's  mass  can  quickly  give  the  minimum  count  time  needed. 

This  problem  will  be  based  on  the  source  spectrum  of  B.  The  source 
spectra  of  B  was  originally  based  on  a  5:1  mass  ratio  of  U:Pu.  From  the 
relative  amounts  of  the  source  spectrum,  a  source  intensity  may  be  found. 
Table  XI  shows  the  finding  of  the  total  intensity  of  Source  B. 


Table  XI.  Total  intensity  determination 
of  source  B 


ENERGY  SOURCE  INTENSITY  EFFECTIVE 


BIN  (KeV) 

SPECTRUM  (%) 

MATERIAL 

(1/g.  sec.) 

INTENSITY  (1/g. sec) 

450 

.1043 

_  239 

Pu 

3.4E4 

3.55E3 

400 

.100 

Pu239 

3.4E4 

3.4E3 

300 

200 

.796 

235 

U  .  6 

4.3E4 

2.85E3 

100 

Pu  ^ 

1.4E5 

1.86E3 

TOTAL:  5.4E4/g.  sec. 

The  minimum  count  time  can  be  determined  by  the  following 
equation. 


Count  Time  =  ‘ 

I^  *  Area 


where 


(1) 


I^  =  Total  intensity  of  source  in  l/g-sec. 

Area  =  Area  of  detector  face 
For  source  B  I^  =  5.4E4  photons/g  sec. 
substituting  into  (1) 

Count  Time  =  ^232-(Range),^  g-  sec. 

Area  ^ 

2 

Figure  35  plots  the  count  time  vs.  range  for  a  detector  area  of  5  cm 


I 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

The  concept  of  utilizing  entire  spectral  features  for  source  iden¬ 
tification  has  been  investigated.  The  study  used  a  simple  correlation 
function  and  employed  two  decision  rules  for  determining  the  source  of 
the  unknown  spectra.  Two  scenarios  were  investigated:  (1)  each  mea¬ 
sured  spectrum  has  a  source  in  the  library,  and  (2)  the  measured  spec¬ 
trum's  source  is  not  a  library  source. 

When  the  decisions  rules  are  applied  to  scenario  (1)  the  following 
conclusions  are  drawn: 

(1)  .  The  counting  of  collided  photons  as  well  as  uncollided  pho¬ 
tons  seems  to  add  structure  to  the  spectrum  which  is  unique  to  the  source. 
The  added  structure  extends  identification  further  than  the  photopeak 
identification. 

(2)  .  The  proposed  techniques  for  identifying  unknown  sources  have 
merit  as  a  possible  alternative  to  photopeak  identification  based  on 
their  success  for  choosing  the  correct  sources. 

(3) .  The  "matching"  of  a  statistically  "poor"  spectrum  still  gives 
results  that  are  correct  in  most  cases  studied. 

When  the  decisions  rules  are  applied  to  scenario  (2)  the  following 
conclusions  are  drawn: 

(1) .  The  small  values  of  the  measured  coefficients  "might"  indi¬ 
cate  that  the  measured  source  is  not  a  library  source.  But  these  low 
values  are  indicated  for  ranges  within  the  photopeak  range. 

(2)  .  Based  on  the  decision  rules  the  non-library  spectra  cannot  be 
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discriminated  as  a  "new"  source. 

Lastly,  a  minimum  count  time  was  determined  for  source  B  as  a  func¬ 
tion  of  range. 

Recommendations 

Based  on  the  assumptions  and  observations  while  working  on  this 
thesis  the  following  recommendations  are  proposed: 

(1)  .  The  modeling  of  a  more  realistic  detector  with  appropriate 
response  functions  should  be  considered.  However,  the  use  of  Morse  for 
photon  transport  is  not  suggested  because  of  its  complexity.  The  only 
reason  Morse  should  be  used  is  to  model  complex  3-D  geometries  which  can 
not  be  done  by  more  conventional  transport  codes. 

(2) .  One  of  the  underlying  assumptions  throughout  this  study  has 
been  that  a  detector  which  allows  both  the  collided  and  uncollided  pho¬ 
tons  produces  a  spectrum  which  can  be  used  for  extending  identification. 
As  mentioned  most  detection  schemes  want  to  reduce  background  (scattered 
photons  included)  and  only  get  peak  structure. 

A  study  of  possible  detector  types  and  designs  which  could  be  used 
should  be  investigated. 

(3) .  Using  the  simple  correlation  equation  on  the  group  structure 
of  the  spectra  was  a  managable  task.  Applying  the  proposed  techniques 
to  real  data  would  involve  a  complicated  process.  This  process  would 
consist  of  setting  up  a  data  acquisition  system  (a  personal  computer 
might  suffice  but  it  would  be  limited  by  storage  capabilities)  and  the 
proper  analysis  routine  (computer  program  to  solve  for  the  coefficients) . 
The  analysis  system  would  involve  finding  the  proper  correlation  tech¬ 
nique  to  handle  the  large  amount  of  data  due  to  the  large  number  of 


channels  In  todays  multichannel  analyzers,  l.e.  finding  coefficients  on 
a  group  to  group  basis  would  be  extremely  costly  in  computer  resources. 
This  work  would  lay  the  foundation  for  further  nuclear  spectral  analysis 
studies. 

(5).  The  entire  analysis  used  a  3  Source  library.  Does  the  iden¬ 
tification  range  increase  as  the  library  becomes  larger?  Do  the  decision 
rules  still  compare  with  one  another  in  their  success  rate?  These  ques¬ 
tions  are  the  basis  for  further  study  for  determining  correlation  limi¬ 
tations  for  the  proposed  decision  rules. 
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Appendix  A.  Photon  Energy  Group  Boundaries 
for  DLC-48/PVC.  (Ref  ) 


Group  Number  Energy  (MeV) 


ii 


I 

r 

n 

i  %  • 

i 


1 

14.0 

2 

12.0 

3 

10.0 

4 

8.0 

5 

7.5 

6 

7.0 

7 

6.5 

8 

6.0 

9 

5.5 

10 

5.0 

11 

4.5 

12 

4.0 

13 

3.5 

14 

3.0 

15 

2.5 

16 

2.0 

17 

1.66 

18 

1.50 

19 

1.33 

20 

1.0 

21 

0.80 

22 

0.70 

23 

0.60 

24 

0.512 

25 

0.510 

26 

0.45 

27 

0.40 

28 

0.30 

29 

0.20 

30 

0.15 

31 

0.10 

32 

0.075 

33 

0.060 

34 

0.045 

35 

0.030 

36 

0.020 

0.010 
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Appendix  B 


A. 


Variance  of  Library  Coefficients, 

2 

The  variance,  a.  „  ,  is  determined  by 
A,B 


2 

AB 

the  following: 


Recall 


But  actual 


Pab 


^1 


+ 


A 

1 


+ 


t. 


B 

1 


where 

=  estimated  means 
=  true  means 
=  random  error 


The  expectation  of  p.„  is 

AB 

j  +  j 
Pab  ■  ®  I  E  1 


now. 


E 


and  E 


=  0 


(1) 


Because  it  is  equally  likely  that  the  different  will  be  positive 
as  well  as  negative,  therefore  the  average  (or  expectation)  will  be  zero. 
(1)  Therefore  reduces  to 


»A.B-  +  (2) 

AR 

if  A=B  then  6^=1 
AR 

if  Aj^B  then  6^  =0 


82 


i.j  i.j 


^AB  • 


M^AB  =^ 


i 


Case  A?^B 


EsA-B-A-B  ^  V*  fn^V  B" 

CiCiCjC.  +  2^ 

i.j  i 

4.  V  a2  ^  a2  b2 


Case  A=B 


i.j  i.j 

+  4  2  C*C^crf  +  '2  of  of 


Therefore 


Case  A^B 


2(*’  •  x  (?)■•;■  •  2  i-iW 


Case  A=B 


2  V  -A^  a2 

'^AA  =  ^  ^1 


The  standard  deviation  is  .  The  standard  deviations  are  actually 
standard  deviations  of  the  means  since  the  a^s  where  of  the  mean.  This 
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is  just  for  the  library. 


B.  Variance  of  Library  Coefficients,  o  „  .  . 

- ^ -  pM,L 

There  are  two  types  of  error  to  consider; 

(1)  The  standard  error  of  the  mean  associated  with  estimating  the 
fluence  values,  i.e.  where  is  the  random  error  of  the 


(2)  The  error  associated  with  the  measured  fluence  C^,  which 
indicates  the  deviation  of  about  C^. 

To  find  o„  ^ 


i 

W-  •  E 

This  reduces  to 


El2  m2 

Oi 
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Therefore 


ML 


■  I  ^  I  ^  I 


Now 


;;M 


C^  must  be  estimated  by  assuming  a  source  L  such  that 


since  no  estimate  of  is  available,  i.e. 


^  {  from  single  1  -  ;;L 
1  y  batch 


while 


/  avg  estimate  after 
i  y  many  batches 

2  2 

In  a  similar  fashion  o^  must  be  estimated.  Since  o^  corresponds 
to  the  variance  of  the  data  and  not  the  mean,  the  relationship  between 
the  assumed  library  source  is 

°i  N 

where  N  =  number  of  batches  used  in  estimating  the  fluence  mean. 

For  the  1000  photon  case 


M 

“l  ■“<’1 


and  for  the  100  case 


m2 

=  10  *  N  *  0^ 
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Since  counting  statistics  says  that  the  estimates  should  vary  as 


; 


I 

. 

I 


B 


B 
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Appendix  C.  Distribution  of  measured  estimates  about  their 
library  estimates. 


LEGEND  FOR  APPENDIX  C 


Bin  (KeV) 

Library  (Counts/Photon) 
Standard  Deviation 
of  the  Mean  (s.d.m.) 

Bin  (KeV) 


Measured  Number  of 
N/A  s.d.  values 

are  from  library 
value 
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H 


8AN6E«41  H 

HEASURED  A- 1000 
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'min« 
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3.ei4E-e06 

8.008E4080 

18 

1.489E-005 

0.88BE40e8 

e.eaeE^aBB 

8. 80084088 

8.8888^800 

0. 08084088 

DEVIATION 
FRON  NEAN 

e.flaBEtaee 

6.723E-aBl 

7.B75E-BBI 

•.24&E^BB1 

&,23BE-BB1 

5.B2BE-BB1 

'.MBEtBBB 

•.1&4E^BB1 

9.33BE-BB1 

-.326EiBBB 

2.179E-BB1 

4.347E-eBl 

6.6S6E-BB1 

1.779EtBBB 

-.363E*BeB 

1.B5BE-BB1 

-.IISE+BBB 

2.562E-BB1 

B.BBBE^BBB 


CGUNTS/PHOTON 
GROUP  LIBRARY  A 
1330 

0.eNE*e00 

e.ePBE^Bae 

1000 

1.914E-002 

4.429E-003 


700 


1.127E-002 

3.a&BE-003 


8.123E-002 

6.779E-00! 


512 


2.880E-002 

5.770E-a03 


510 


3.555E-802 

3.301E-0a3 


3.514E-0 


450 


'S.487E?iB3 


400 


3.432E-002 

7.147E-003 


300 


8.816E-002 

l.!01^002 


200 


1.234E'001 

2.135E-002 


150 


9.955E-002 

2.160E-002 


100 


1.635E-001 

2.622E-e02 


1.172E-001 

3.273E-002 


60 


7.679E-002 

2.0108-002 


45 


8.765E-002 

2.939E-002 


30 


4.260E'002 
1.933E-002 


20 


2.912E-803 

2.051E-003 


10 


1,739E-B06 

S.440E-006 


g.0B0E+000 

0.000E^000 


RAN6EM66  tt 
HEASUREO  A- 1000 


0.000E>000 

8.000Et880 


1.747E-002 
0. 0008^000 


1.541E-0., 
0.0a0E*i00 


8.1S4E-002 

0.000E+000 


1.637E-002 

0.00BE*080 


3.283E-002 

0.000Et008 


3.748E-B02 


a.asBE^iaa 


3.121E-002 

0.a00E^0aa 


9.e56E-a02 

0.000E^00a 


1.199E-a01 

0.000E>000 


1.335E 

i.r 


0,000MB0 


1.444E-001 

0.00BE*808 


9.454E-I02 

0.000^000 


5.357E-B02 

0.000E^000 


1.196E-001 

0.000E^000 


3.558E-802 

0.000Ef0B0 


4.286E-003 

0.800E^008 


0.0008^000 

0.000Et000 


0.000Et00|. 


0.000E*800 


DEVIATION 
FROH  HEAN 


0.000E+000 

3.762E-001 


-.I07E+001 


-.465E-001 

2.1S3Ei00O 

8.2SIE-001 

-.4118^000 

4.3S2E-001 

-.6938^000 

1.600E-001 

-.1S7E*001 

7.269E-001 

6.932E-001 

1.1558^000 

-.ia8E+001 

3.676E-001 

-.669E+000 

3.197E-801 


0.0008*000 


I. 


I 


CQUNTS/PHOTON 

RAN6E>24S  K 

DEVIATION 

GROUP 

1330 

LIBRARY  A 

HEASURED  A-1000 

FROn  HEAN 

0.000E+000 

0, 0008^000 

0.00084000 

1000 

0.000E*000 

0,00084000 

1.074E-002 

1.161E-002 

-.40284000 

G00 

2.1S2E-003 

0.00084000 

6. 137E-003 

1.101E-002 

-.24084001 

700 

2.024E-003 

0,000E4000 

4.0B7E-002 

3.114E-002 

1.68784000 

600 

5.292E-003 

0,000E40B0 

2.041E-002 

2.263E-002 

-.48584000 

512 

4.566E-003 

0,0B0E4000 

1.438E-002 

1.353E-002 

4.2968-001 

510 

1.957E-003 

8. 00084080 

2.352E-B02 

2.0868-002 

5.4628-001 

450 

4.a76E-003 

0,00084000 

2.371E-002 

2.1388-002 

5.0728-001 

400 

4.599E-003 

0,00084000 

6.089E-002 

5.8438-002 

2.0068-001 

300 

200 

1.226E-002 

0.00084000 

9.273E*002 

1.944E<802 

'»« 

6.7298-001 

6.482E-002 

5.6028-002 

8.6828-001 

150 

1.0148-002 

0.00084000 

1.292E-B01 

1.3888-001 

-.41784000 

1  sa 

2.283E-0B2 

0.00084000 

IIBQ 

9.B84E-002 

8.6528-002 

2.2138-001 

75 

1.954E-002 

0.00084000 

7.496E-002 

4.9398-002 

6.2128-001 

60 

4.116E-002 

0.00084000 

8.839E-002 

6.2208-002 

8.5918-001 

45 

3.049E-002 

0.00084000 

4.720E-002 

2.0338-002 

1.28384000 

30 

1.471E-002 

0.00084000 

4.174E-003 

5.6728-004 

7.0468-001 

20 

5,ll9E-003 

0.00084000 

1.234E-B06 

0.00084000 

2.8988-001 

10 

4.257E-006 

0.00084000 

0.000E*000 

0.00084000 

0.00084000 

0.000E*000 

0.00084000 

94 


I 


COUNTS/PHOTOM 
GROUP  LIBRARY  A 
1338 

8.ee8E*B8e 

e.eeeE^See 

1888 

5.488E-ee3 

1.187E-B83 

888 

3.5838-883 
1.1888-883 


788 


688 


512 


518 


458 


388 


288 


158 


75 


be 


45 


38 


28 


18 


1.985E-882 

2.275E-883 

i.l32E-eB2 

2.533E-ee3 

5.4888- 883 
1.3658-883 

1.4858-882 

3.2748-883 

1.3468-882 

3.7878-883 

3. 9598-882 
5.2978-983 

6.3888-882 

3.1188-882 

4.2488-882 

9.8918-883 

8.3368-812 

1.6998-882 

6.4558-082 
1.3538-082 

6.1478-082 

2.3818-882 

6.4888- 882 
2.0578-082 

3. 1868-802 
1.1928-002 

2.3768-003 
2. 4678-003 

5.4738-006 
1.8118-085 

8.8808^800 

0.8088^080 


RAN6E-333  N 

HEASURED  A-IBBS 

0.ee0E>e0e 

e.ee0Et880 

4.8638-003 

0.0008+008 

4.7208-003 

0.0808+000 

2.1988-002 

0.0808+080 

1.4948-002 

0.0008+008 

4.1318-003 

8.8808+880 

1.3748-002 

8.0808+008 

2.8368-002 

8.0008+888 

5.5788-002 

8.0888+000 

5.1088-802 

8.8808+000 

3.2638-082 

8.8008+088 

9.0768-002 

0.0808+888 

5.6078-002 

8.0888+000 

9.7808-002 

0.0088+008 

1.2198-001 

8.0088+000 

4.6008-802 

8.0008+000 

1.0068-883 

8.8808+888 

8.8808+800 

8.0008+808 

8.8008+088 


DEVIATION 
FROM  KEAN 


0.0008+000 
5.2638-081 
-.  1838+001 
-.9388+800 
-.1438+001 
9.8848-081 
3.3948-001 
-.1868+801 
-.3048+081 
4.1188-001 
9.9628-001 
-.4358+088 
5.6768-001 
-.1548+001 
-.2888+801 
-.1188+081 
5.5548-001 
3.0238-001 
8.0888+888 
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COUNTS/PHQTON 

RANGE<390  K 

DEVIATION 

GROUP  LIBRARY  A 
1330 

HEASURED  A-1000 

FROn  HEAN 

3.0458-103 

B.BBBE^BBB 

0.08IE400B 

0.00084008 

10B0 

3.4B7E-0B3 

7. 0748-003 

1.2SaE-8B3 

0.00084000 

-.29184001 

808 

2.389E-0a3 

1.873E-B03 

^'a!lS8!4i08 

-.87084000 

700 

t.l72E-8B2 

5.4648-003 

2.8258-003 

0.00084000 

2.21584000 

600 

8.181E-003 

7.4188-003 

2.723E-003 

0.00084000 

2.8018-001 

512 

2.726E-003 

0.00084000 

6.445E-004 

0.00084000 

4.23084000 

510 

9.3828-003 

4,6428-003 

2.233E-003 

0.00084000 

2.34784000 

450 

1.0358-002 

3.3848-003 

3.952E-003 

0.00084000 

1.76284088 

400 

2.S84E-802 

2.3818-002 

2.9528-001 

6.896E-003 

0.00084000 

300 

3.8618-082 

3.6328-002 

2.3078-001 

9.940E-003 

0.00084000 

200 

3.106E-082 

1.043E-002 

3.3148-002 

0.00084000 

-.19984000 

150 

6. 159E-002 

4.6198-002 

2.1488-002 

0.00084000 

7.170E-0B1 

100 

4.9418-002 

6.1188-802 

2.3038-002 

0.000^000 

-.55384000 

75 

4.0078-002 

5.5928-002 

1.183E-002 

0.80084000 

-.13484001 

60 

5.2028-002 

3. SI  38-082 

1.602E-002 

0.00084008 

8.6758-001 

45 

2.239E-002 

2.3698-002 

5,7728-003 

0.00084000 

-.22584008 

30 

2.4408-003 

3.9848-084 

2.2248-003 

0.00084000 

6.9308-001 

20 

9.. .398-006 

0.00084000 

2.4668-005 

0.00084000 

3.4018-001 

10 

0.3008t000 

0.00084000 

0.0008t000 

0.00084000 

0.00084000 

e.eeeE^eze 


COUNTS/PHOTON 
GROUP  LiBRftRV  B 

^^’'SlftSUREO  B-1000 

B8VlftT10M 
FROn  HERN 

0.000Et005„ 

0.0008+000 

*-!We800 

0.0008+000 

‘-rBau 

0,0008+000 

0.0008+000 

0.0008+000 

800 

0.0008+000  . 
0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

790 

0.8008+000  . 
0.0008+000 

0.8008+000 

0.0008+000 

0.0008+000 

400 

0,0008+080,. 

0.0008+000 

0.0008+000  , 
0.0008+000 

0.0008+000 

0.0008+000  , 
0.0008+000 

8.0008+000 

0.0008+000 

0.0008+000 

5ia 


m 


200 


100 


75 


te 


30 


20 


10 


0.00eE>i00 


7.‘)«E-002 

l.003E-0a2 

9.4i3£-002 

1.057E-002 

2.445E-002  , 
1.124E-002 

i,406E:0ll  , 
l.(i77E-002 


1.759E-001 

4.177E-002 

3.174E-0I2 

l.502E'002 

4.879E-003 

ti.7a8E'003 

1.660E-003 

4.435E-0a3 

0,000E»000 

0.000E^000 

0.0000^000 

0.0008^000 

0.000E*000 
0.0008^000 

.0008^000 
0.0008^000 


9.359E-002 

0.0008^000 

1.024E-001 

0.0008^000 

1.4fa3£-002 

0.0008^000 

4.5778-001 


'’H00S000 

3.1658-002 


1.3518-002 
0. 0008^000 

1.2798-002 

0.0008^000 


0.0008^000 

-.1418^001 

-.7548+000 

8.7358-001 

1.7578-001 

4.1918-001 

-.7488-002 

-.1278+001 

-.2508+001 


0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

0.0008+000 

97 


CBUNIS/PHOTON 

GROUP  LIBRARY  C 
13J0 

B.eeeE+aeB 
B.eaeEteae 


2.2esE-eBi 

l.A5eE-0e2 

i.i2iE-0ei 

9.743E-003 

3.0nE-002 

5.3&2E-003 

8.795E-0B3 

4.728E-003 

4.417E-004 

1.042E-803 

A.838E-003 
4.270E-003 

5.9n£-003 

3.7B7E-003 

1.29&E-0e2 

8.871E-003 

2.UeE-002 

3.983E-003 

I.lli8£-001 

1.817E-002 

1.536E-001 
2.171E-082 

4.689E-002 
1.228E-002 

8.653E-083 
7.208E-003 

2,327E-003 

5.245E-003 

3.638E-004 
1.6258-003 


700 


600 


512 


510 


450 


400 


300 


200 


:50 


100 


75 


60 


45 


30 


0. 0008*000 


20 


RAR6E*2e.  N 

HEASURED  C-1000 

0.000E*00i 

0.000E*000 

2.261E-001 

0.000E*000 

1.194E-001 

0.000E*800 

3.605E-002 

0.000E*000 

1.630E-0B2 

0.000E*000 

2.129E-083 

0.000E*000 

8.227E-003 

i.000E*880 

1.374E-803 

0,0008*080 

1.941E-082 

0.800E*000 

2.651E-002 


10 


0.0008*000 

0.0008*000 

0.0008*000 


4.354E-001 

0.0008*000 

1.8608-001 

0.0008*000 

4.0188-002 

0.0808*000 

9.8368-003 

0.0008*000 

0.0088*000 

0.0008*880 

0.0088*080 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 


DEVIATION 
FROn  HEAN 

0.0808*000 

-.3888*080 

-.7528*000 

-.1068*001 

-.1588*001 

-.1628*001 

-.3258*000 

1.1988*000 

-.7998*000 

-.6028*000 

-.1028+001 

-.1498+001 

5.4608-001 

-.1718*000 

4.4368-001 

2.2398-001 

0.0008*000 

0.0008*000 

0.0008*000 


98 


B  %  * 


CQUNTS/PHQTON 
GROUP  LIBRARY  B 

RAN8E=a3  K 

H8ASURED  8-1000 

DEVIATION 
FROM  H8AN 

0.0008*000 

0.0008*000 

0.0008*800 

0,0008*000 

0.8008*000 

1000 

0.0008*000 

0.0008*000 

0.0008*000  ^ 
0,0008*000 

0,8008*000 

800 

1.6518-005 

9.0448-005 

0.0008*000 

0.0888*000 

1.8268-001 

700 

0.0008*000 

0.0008*000 

'I'HitnH 

0,0008*000 

600 

0.0008*000 

0.0008*000 

^'0J00l*0B0 

0.0008*000 

0.0008*000 

0.009Ei8««  . 

a  Baac^BnB 

510 


450 


400 


300 


200 


150 


100 


75 


!>0 


30 


0,000E>000 

0.000E*000 

4.l<!9E-002  , 
6.078E-003 

6,44lE-002  , 
7.874E-003 

4.505E'002 

1.179E-002 


3.559E-801 

l,779E-002 

3.489E'00i 

2.4UE-002 

1.626E'001 
3.312E-802 

a.667E'002 

1.879E-002 

5.630E'002 

2.057E-002 

K718E-002 

1,172E'002 

3.792E-004 

1.422E-003 


0. 

0.000Em0 

0. 0008^000 

4.l5flE:M: 


5.040E'002 
0. 0008*000 

*’0!800!*i80 

^'M00E*000 

3.7B8E-001 

0,0008*000 

1.9128-011 

0.0008*000 

8.2098-002 

0. 


10 


'0.0008*000 

0.0008*000.. 

0.0008*000 


5.7498-002 

0.0008*000 

2.1718-002 

0.0008*000 

9.8958-005 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 


0008*000 

7458-002 

7798*000 

6,0348-002 

5.1388-001 

-.1578*001 

-.8638*000 

2.4418-001 

-.5778-001 

-.3868*000 

5.4898-001 

0.0008*000 

0.0008*000 
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COUNTS/PHOTON 

RANGE: 124. 5  H 

DEVIATION 

GROUP  LIBRARY  8 

MEASURED  8-1888 

FROn  MEAN 

1338 

e.eeBEveee 

8.8888*888 

a.eeeEtaee 

0.0088*888 

e.800E*000 

1888 

B.eeBE^eae 

8.8888*008 

8.8eeE«eee 

0.0808*800 

e.008E*00e 

388 

e.e8eE>ee8 

8.8888*088 

B.eeeE^eea 

0.8888*888 

8.0888*888 

788 

e.e88Et88e 

0.8088*888 

8.8eeE>88e 

8.0888*888 

8.8888*808 

388 

8.8eBEtee8 

8.8008*080 

e.888E^888 

8.8888*888 

0.0808*088 

j  1 2 

8.8888^888 

8.8888*808 

8. 8888^388 

8.8888*888 

8.8888*808 

;i8 

8.8888*888 

8.8808*000 

8.8888*888 

8.8808*888 

0.8888*888 

458 

2.78')E-882 

2.3028-002 

3.6568-883 

8.8008*888 

-.255E+000 

488 

4.5658-882 

5.0928-002 

8.8628-883 

8.8888*888 

-.6538*008 

388 

4.7848-882 

4.4138-882 

1.5498-882 

8.8808*008 

2.363E-801 

288 

2.2828-881 

2.3838-001 

1.5898-802 

8.8888*888 

-.6318*800 

158 

3.1648-081 

2.3518-001 

4.5808-802 

3.8888*888 

6.933E-001 

138 

1.3628-881 

1.7128-001 

8.8888*808 

4.8228-082 

3.1058-801 

1.1678-801 

1.3318-801 

2.8198-002 

8.8008*888 

3.3468-001 

68 

1.1448-001 

1.8468-881 

4.4518-802 

8.8088*088 

2.196E-001 

45 

^ 4468-082 

2.7728-802 

1.5278-802 

8.8088*888 

7.6908-801 

30 

3. 1588-083 

1,9098-004 

3.2078-803 

8.8088*888 

9.25*8-001 

28 

1.3628-005 

1.5518-007 

5.1318-805 

8.8888*888 

2,6648-031 

18 

0.8008*388 

0.8088*808 

0.3888*888 

3.8888*888 

8.0808*808 

COUNTS/PHOTON 
SROUP  LIBRARY  B 
i3JB 

B.BME^Ne 

e.BBBE^BBI 

IBBB 

e.enE^eii 

B.eSeE^Bai 

OAA 

RAN6E>1&6  H 

HEASUREO  B-ieea 

DEVIATION 

FRDN  HEAN 

a.BaaE^eBe 

a.eaaE^eia 

‘SIm 

e.aaeE+eee 

a.eaeEma 

QvV 

B.BBIE^MB 

B.eeeE^flea 

7ee 

e.fle«E>eae 

e.flaaE^BBB 

6es 

7.a29E-«a5 

3.8SeE-Be4 

512 

e.aaaE^aea 

e.eeeE^eee 

5ia 

a.eaaE^aea 

B.eaeE^eea 

45a 

i.62iE-ea2 

2.523E-8a3 

4Ba 

3.568E-ae2 

7.e84E-aa3 

3ee 

4.356E-aa2 

i.ai5£-a82 

208 

i.492E-eei 

1.392E-002 

isa 

2.427E-001 

2.7l7E*a02 

\  Afl 

B.eeeE^eai 

e.eeaE^eaa 

e.aeeEieee 

e.aaeE^eaa 

e.eaaE^aae 

a.aBBEma 

e.eaeE+aaa 

a.eaaE^aeB 

i.026E-eei 

a.BeeE^eae 

a.eeaE^BBa 

B.eaeE^BBa 

a.eaeE^aaa 

a.eeaEfeea 

0.a88£fB00 

'‘MSeliiaa 

-.&S9E^88a 

3.80i£-ae2 

e.00eE*0a8 

-.329E+000 

4.e48E-8B2 

B.00aE>000 

-.4&2E*eee 

- - - 

1.515E-0ai 

B.BBiEiaea 

U5Ete00 

-  i 

2.927E-001 

a.BBBE^aaa 

•.184E^881 

1  00 

1,757E-001 

3.365E-a02 

75 

i.mE-eai 

2.968E-002 

60 

1.259E-001 

3.78iE-a82 

45 

5.949E-002 

l.866E-ee2 

30 

4.316E-ea3 

2.371E-003 

20 

2.702E-006 

1.2138-005 

10 

0.00aE^000 

0,eaeE*e00 

1.751E-001 

a.eaaE^eea 

1.785E-002 

•  1 

1.235E-I01 

e.aSaE^a 

1.369E-881 

‘1 

l,043E-e81 

e.000E>a00 

5,7I4E-001 

^■0!?l8!S^00 

5.037E-801 

1.294E-003 

0.000E^ee0 

1.274Ete08 

—  .  • 

8.402E-007 

a.aaeE^aaa 

1.535E-001 

o.eeBE^eaa 

a.aeeE^aaa 

0.000E*88e 
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CDUNTS/PHOTON 

RAN6E-245  K 

DEVIATION 

GROUP  LIBRARY  B 
1338 

i.7eiE-ees 

9.3ieE-ies 

iflee 

B.eeeE^eee 

MEASURED  0-1000 

FROM  MEAN 

e,00eE*00B 

0. 0000*000 

1.824E-001 

0.000E*000 

0.00BE*000 

0.000E*B00 

0vO 

B.eeeE^BBe 

e.eeeE^BBe 

0.000E*B00 

0.0000*000 

8.C00E*000 

700 

0.000E+00e 

0.e0eE*000 

ma 

0.000E*000 

8.000E*000 

0.0000*000 

olfliu 

0.0eeE*0B0 

0.e00E>0e0 

0.0000*000 

B.B00E*000 

0.000E*0B8 

f  12 

0.000E*eBe 

0.000E>0e0 

510 

0.0eeE>e0e 

0.000Ei000 

0.000E*08B 

0.0000*000 

0.0000*000 

0.00BE*000 

0.0000*000 

0.000E*000 

450 

6.aeeE-0a3 

1.804E-003 

4.603E-003 

0.00BE*000 

1.539E-001 

400 

1.753E-002 

3.246E-003 

1.295E-002 

0.088E*000 

1.4130*000 

2.650E-002 

6.524E-003 

200 

4,748£*002 

1.1  BSE -002 

3.184E-002, 

0.0000*000 

-.618E*B00 

6.864E-002 

0.0000*000 

7.529E-002 

150 

1.379E-001 

1.7BiE-002 

1.557E-001 

0.000E*B00 

-.99SE*000 

ilSil 

l.lllE-001 

l.a42E-002 

75 

9,930E-002 

2.100E-002 

b0 

1.124E-001 

2.495E-002 

45 

5,875£-002 

1.447E-002 

1.139E-001 

0.0000*000 

-.151E*000 

6.S2iE-002 

0.0000*000 

1.478E*000 

1,220E-001 

8.000E*000 

-.4170*000 

8. 8300-002 

0.000E*000 

-.Z04£+00I 

3.274E-003 

2.193E-003 

20 

1.766E-005 

4.446E-005 

10 

0.000E^000 

0.00eE>000 

4,5420-003 

0.0000*000 

-.578E*000 

0.0000*000 

0.0000*000 

2.740E-001 

0.0000*000 

0.0000*000 

0.0000*000 
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c# 


COUNTS/PHQTON 

RANGE:333  H 

DEVIATION 

GROUP  LIBRARY  8 

HEASURED  0-1000 

FROH  HEAN 

1330 

0.e00Eteee 

0.0eaE>i0B 

4  (lAfl 

0.080E4000 

0.00aE400B 

0,0008*000 

1  O0V 

0.000E+e08 

0.000Eteee 

flAA 

0.0008*000 

0.000E*000 

0,0008*000 

OUv 

0.0B0E*eBe 

B.eeeE^eoa 

*TAA 

0.000E*000 

0,000E*000 

0,8088*000 

700 

0.eeeE>0a0 

B.BaBE^aae 

1  Art 

0.0008*000 

0.0008*000 

0,0808*000 

000 

0.80BE^0e0 

B.aaeEtaee 

0.0008*000 

0,0008*000 

0,0008*000 

512 

e.eeaE^eai 

e.eaaE^eaa 

an 

0.0008*000 

0.0008*000 

0,0008*000 

510 

0.e08E*eii 

e.eaaE^aea 

450 

:,635E-003 

1,139E-B03 

0.0008*000 

0,0008*000 

0,0008*000 

2.5798-003 

0,0008*000 

4,8618-002 

400 

7.894E-003 

1.909E-003 

9.7778-003 

0.0008*000 

-.9868*000 

300 

1.487E-002 

5.2488-003 

200 

2.836E-002 

6.912E-003 

150 

6.B84E-002 

1.449E-002 

100 

5.934E-002 

1.883E-002 

1.4868-002 

0.0008*000 

1.7968-003 

4.0038-002 

0.0008*000 

-.1688*001 

6.5188-002  . 
0.0008*000 

2.5428-001 

5,3268-002 

0.0008*000 

3.2288-001 

75 

4.8708-002 

1.225E-002 

/  A 

5.8178-002  , 
0.0008*000 

-.7738*000 

00 

6.996E-002 

2.029E-002 

45 

4.081E-002 

2,il7E-0a2 

6.2038-002 

0.0008*000 

3.9088-001 

3.5078-002 

0.0008*000 

2.1918-001 

30 

3.007E-003 

2.540E-003 

1.4808-003 

0.0008*000 

6.0098-001 

*'0 

6.250E-004 

2.233E-005 

0.0008*000 

0.0008*000 

2.7988-001 

10 

0.008E^000 

0.0008^000 

0.0008*000 

0,0008*000 

0.0008*000 
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CQUNTS/PHOTON 

range-390  n 

DEVIATION 

GROUP 

'  LIBRARY  B 

MEASURED  B-1800 

FROM  MEAN 

1338 

8.e88E68ee 

e.eeoEieee 

e.8e8E«eee 

e.aeaE^oea 

0.0008*008 

1082 

e.eeeEtaee 

0.0e0E60ee 

9.98flEteee 

0.e00E*8ee 

0.00084000 

see 

8.eeeE+aae 

e.eaeE^eee 

e.eeeE^eee 

0.800E*000 

0.0008*000 

788 

e.eaeE^aee 

0.0008^000 

a.eeeE^aea 

0.0808^800 

0.0008*000 

^88 

e.eeeEieee 

8.0088*008 

511 

8.eeeE>888 

8.8008*800 

8.8008*000 

e.eeeE^eee 

8.8008*080 

a.aaaE^aea 

0.0088*000 

0.0008*000 

5!fl 

a.eaaE^eee 

0.0008*000 

a.eeeE^eee 

8.8008*880 

0.0008*000 

450 

1.367E'aa3 

1.0318-003 

5.974E-e84 

0.8008*800 

5.6368-001 

408 

4,655E-883 

4.203E-003 

1.384E-803 

0.0008*000 

2,2658-001 

328 

8.349E-083 

9.6208-003 

3.524E-ee3 

0.0008*808 

-.4178*000 

288 

1.436E-082 

1.0986-002 

150 

4.205E-e03 

0.0008*808 

8.0528-801 

j.aisE-aa: 

O.736E-002 

1.8888-882 

0.8008*808 

-.2688*001 

188 

3.382E-002 

4,3498-082 

7.9848-883 

0.0808*800 

-.183E*001 

75 

3. 121 £-802 

3.295E-002 

7.35lE-0e3 

0.0088*808 

-.236E*000 

68 

4.3448-882 

4.65BE-082 

1.2878-082 

8.8808*080 

-.399E*:O0 

2.5218-802 

2.4928-002 

3.718E-003 

8.0008*800 

3.3248-002 

TJI 

!.:'41E-B03 

9.4738-004 

28 

1.7648-003 

0.0008*080 

4.5008-001 

1.3=38-006 

0.0008*000 

1  n 

6.5478-006 

0.0008*088 

2.1278-801 

1 V 

8.0008^000 

0.0808*800 

e.eooE^eea 

0.0008*000 

0.0008*000 
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i 

I 


COUNTS/PHOTON 
GROUP  library  B 
133B 

a.eBBE^Bee 

e.flBeEfBBe 

leee 

B.BBBE^BtB 

B.BaBE^eea 


aea 


512 


j.eeaE^eea 

a.eeaE^aee 

a.eeaE+eee 

a.eaaE^eee 

e.aeaE^aaa 

e.eaeE^aae 

e.eaaE^eae 


sie 


450 


400 


m 


200 


150 


100 


75 


60 


45 


30 


20 


10 


0.000E*00B 


6.4&0E-602 

8,0S9E-B03 

7.944E'002 

1.104E-a02 

4.009E-082 

1.351E-002 

5,425E-001 

2.573E-002 

2.726E-001 

3.312E-002 

7,534E-002 

2.373E-002 

3.223E-B02 

2.258E-002 

1.045E-BB2 

1.145E-802 

2.06lE-0a3 
3.935E-003 

7.984E-005 
3. 3368-004 

3.091E-006 

1.693E-005 

0.e00E>000 

0.000E^000 


RANGEMl  K 
IIEASURED  B-1000 

0.0BBE*000 

0.0008^000 

3.0e0E>e00 

0.0008^000 

0.000E>000 

0.0008^000 

0.0008^000 


0.00084000 

0.00084000 

0.00084000 

0.00084000 

0.00084000 

0,00084000 

5.7308-002 

0.00084000 

6.7898-002 

0.00084000 

4.2998-002 

0.00084000 

3.6358-001 

0.00084000 

2. 8158-001 
0.00084000 

7.2778-002 

0.00084000 

2.256E-002 

0.00084000 

5.5568-003 

0.00084000 

0.00084000 

0.00084000 

0.00084000 

0.00084000 

0.00084000 

0.00084000 

0.00084000 

0.00084000 


DEVIATION 
FROM  MEAN 

0.00084000 

0.00084000 

0.00084000 

0,00084000 

0.00084000 

0.00084000 

0.00084000 

9.0248-001 

1.04884000 

-.21484000 

-.31784000 

-.26984000 

1.BB2E-001 

4,2788-001 

4.2708-001 

5.2378-001 

2.3938-001 

1,S26E-001 

0.00084000 
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R 


i  if 


COUNTS/PHDTON 
GROUP  LIBRARY  C 
1333 

0.08BE*80B 

0.83eE*88R 

i  snA 

RANBE*41  «  ,  , 
MEASURED  C-1B80 

DEVIATION 

FROM  MEAN 

- 

8.BBBE*BBB 

B.BBBE4BBB 

B.BBBE*BB3 

1 006 

1.832E-001 

1.374E-002 

1.B71E-BB1 

3.00BE*BBB 

-.2a8E*BBB 

366 

9.a36E-832 

3.S60E-083 

AA 

1.003E-BB1  . 
0.flBB£*B00 

-.112E*B01 

- 

706 

3.375E-002 

6.551E-003 

4.4iBE-082 

B.000E*BBB 

-.165E*001 

830 

1.812E-002 

i.l83E-003 

1.896E-0B2 

3.30BE*B00 

-,136E+aBB 

512 

2.563E-804 

1.B33E-0B3 

-,245E*B01 

6.415E-004 

0.003E*030 

510 

l.a35E-002 

5.i29E-003 

1.127E-0B2 

3.000E*0Ba 

-,184E*30B 

450 

8.699E-803 

5.371E-003 

A  .1(1 

7.5S5E-003 

B.0BBE*30B 

2,075E-83l 

400 

r.939E-002 

4.299E-002 

-.123E*00l 

1.129E-002 

0.BBBE*B00 

330 

3.158E-8B2 

1.339E-002 

2.208E-0B2 

B.30BE*030 

7.096E-001 

3.453E-081 

2.111E-002 

3.212E-001 

B.B0BE*BBB 

1,142E*00B 

1 50 

2.209£-a01 

2.786E-002 

2.089E'001 

3.B0BE*BBa 

4.305E-0B1 

!00 

7.805E-002 

2.735E-002 

fc.323E-002  . 
3.000E*B33 

5,418E-B01 

et! 


45 


13 


1 3 


:.?80E-B32 

1.3388-302 
1.010E-002 

3,  <1348-383 
7.066E-i03 

i.336E-035 

:.264£-00< 

e.33eEt8e3 

3.3838*808 

3.0e0E*088 

0.030E*0e0 


..954E-082 


2.343.E-032 

0.000E*380 

!.328£-002 

3.383E*380 

3.330Et008 


3.033Et088 

2. e03E*008 

3.e03E*000 

3. e00E*88e 


2.348E-002 

-.102E+001 

-.139E+801 

2.461E-001 

0.303E*33e 

8.803£*038 
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P 


I 


I 

I 


i 


i 


I 


I 


COUNTS/PHOTON 
3ROUP  LIBRARY  C 
1338 

e.easE^eae 
e.88eE^<fl8 


!,371E-B01 

i.iieE-8e2 


888 

788 

m 

512 

518 

458 

488 

388 

288 

158 


o 


68 


45 


38 


20 


18 


7.732E-082 

8.754E-803 

3.231E-002 

5.545E-a83 

2.094E-002 

4.683E-8e3 

2.727E-a84 

6.2a5£-aa4 

1.624E-002 

7.641E-003 

1.457E-002 

6.220E-003 

3.R6‘)£-8a2 

1.364E-002 

5.082E-082 

1.458E-002 

2.317E-0B1 

1.68aE-002 

2.645E-081 

3.1B3E-002 

1.394E*00l 

2.097E-002 

7.1BBE-B02 

1.969E-002 

5.i63E-002 

1.855E-002 

1.548E-002 

1.220E-802 

9.667E-004 

l.b95E-083 

5.104E-805 

2.795E'884 

8.880E*00e 

e.88eE*080 


RAHBE=03.  M 

HEASURED  0-1888 

DEVIATION 
FROn  MEAN 

e.008E«0B8.. 

e.B0eE*000 

a.aaaEt-Boa 

1,477E-001 

e.000E^0a0 

-.953E^000 

8,075E-a82 

8.00aE^e0a 

-.392E+000 

3.6626-002 

e.0aaE*aB0 

-,776E+e0a 

2.784E-002 

8.000E^000 

-,149E+001 

e.B00E«000 

a.000E^e00 

4.338E-081 

2.176E-002 

e,80aE*0a0 

-.719E+000 

1.29e£-802 

a,00aE*0a0 

2.682E-001 

4.193E-002 

0.e00E*0a8 

-.164E+000 

4.199E-002 

a.8e0E^088 

6.058E-001 

2.323E-001 

8.000E^aB8 

-.39eE-001 

2.68eE-e01 

e.0e0E^000 

1.421E-001 

l,529E-001 

a.000E^0e8 

-.645Ei000 

7.605E-002 

e.eaaE^eea 

-.256Et00a 

3.291E-002 

a.eeaE^aaa 

1.0Z0Ei000 

3.753E-002 

a.e0BE>88e 

-.100E+001 

e.0a8E^88a 

a.e00Et8ae 

5.736E-081 

B.000E^B00 

e.eaeE^aaa 

1.B26E-081 

c.aeeE^aea 

e.8e8E«888 

8. 8880 ^888 

» 


» 
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CQiJNTS/PHQTON 
8RQUF  LlBRftRV  C 
133R 


RANGEM24.S  «  , 

MEASURED  C-1000 


9.933E-B02 

7,370E-003 

5.596E-002 

8.879E-003 

3.078E-0B2 

7.579E-O03 

1.3ab£-092 

i.2a7E-003 

5.999E-804 
9, 7378-004 

'  1.778E-002 
b.l96E-003 

'  1.664E-002 
6.094E-003 

*  4.167E-002 

1.174E-002 

*  5.3O5E-002. 

1.400E-f02 

*  1.626E-001 

I.350E-002 

0 

:'.4B0E-001 

3.281E-002 

« 

!.499t-00! 

2.b63E-002 

^  1.008E-0I1 
2.601E-002 

i 

9.415E-002 

■..4'j3E'B02 

!.37eE-002 

1.3I:E-002 

2.395E-003 
2, 160E'00j 

1.401E-007 

4.565E-B07 

[0 

0.000E+000 

0.e00E‘000 


0.a00E^000 

9.229E-002 

0.009E*00a 

6.102E-002 
0. 0008*000 

2.533E-002 

0.000E*000 

'i.3A9E-002 

0.000E*000 

2.128E-003 

0.0008*000 

i.383E-032 

0.0008*000 

l.fc22E-002 

0,0008*000 

4.0S8E-002 

0.0008*000 

S.8bb£-902 

0.0008*000 

1.940E-00i 

0.0008*000 

2.7598-001 

0.0008*000 

!, 7118-001 
0.000E*000 

1.1468-001 

0.0008*000 

5.26lE-00r 

0.0008*000 

4.4348-002 

0.0008*000 

5.l0b£-003 

0.0008*000 

1.5178-007 

0.0008*000 

0.0008*000 

0.0008*000 


DEVIATICN 
FROM  HEAH 


0.0008*000 


9.5538-001 


-.5698*000 


7.1888-001 


6.0138-002 


-.1578*20! 


6.396E-001 


6.9218-002 


6.7068-002 


-.4008*000 


.1698*901 


-.8708*000 


-.7958*000 


-.5308+000 


4.5868-002 


-.4208*000 


-.1248*001 


-.255E-001 


0.0008*000 


COUMTS/PHOTON 

RAN6EM66  H 

DEVIATION 

6ROUP 

133B 

iBee 

LIBRARY  C 

HEASURED  0-1008 

FROn  MEAN 

0.eNE«0ei 

e.8a0Efi00 

a.aaaE^aaa 

e.aeeEteea 

0.8008+008 

7.19BE-B02 

7.e69E-eB2 

a.aeaEvaaa 

1.668E-B01 

m 

7.289E-003 

4.356E-ei2 

4.696E-aB2 

-.6298+008 

708 

S.4aiE-003 

a.aeBE*000 

2.719E'aa2 

3.589E-a82 

-.130E+001 

X.  AS 

6.696E-0a3 

a.aaaE^aae 

oBi 

2.236E-0a2 

1.77BE-002 

7.6978-801 

312 

510 

6.assE-0a3 

B.B0BE*BB0 

7.361E-ea4 

1.485E-003 

0.BBBE+aBB 

0.8088^088 

4.9568-081 

1.666E-aa2 

2.527E-882 

-.1218+001 

450 

7.107E-8a3 

8.0088+888 

2.03iE-aa2 

3.719E-003 

3.1778-802 

8.0888+880 

-.2808+081 

ABB 

4,i7aE-aa2 

4.527E-002 

-.3098+080 

l.l54E*a02 

8.0808+808 

*00 

5.99lE-0a2 

3.909E-002 

1.2178+000 

200 

1.711E-002 

8.880E+008 

1.241E-0ai 

1.1988-881 

4.0178-001 

( ca 

1.29eE-002 

0.0088+008 

I  wB 

1.947E-081 

1.6828-881 

9.8938-001 

100 

2.672E-a02 

8.e08E+080 

l,476£-0ai 

1.5618-881 

-.2698+000 

75 

3.141E-802 

8.ee8£+008 

1,029E-001 

7.3088-882 

1.30BE+000 

60 

2.279E-002 

0.0008+080 

i.0iBE-0ei 

1. 2348-001 

-.6648+000 

<5 

2.494E-0a2 

0.0088+008 

4,597E-B02 

5.7458-082 

-.612E+00e 

30 

l.fl74E-002 

0.0008+080 

4.412E-e83 

1.431E-083 

6.1508-001 

20 

4.847E-003 

0.0088+000 

4.e45E-e06 

0.0088+880 

3.102E-B01 

1.304E-e0S 

8.0088+000 

IB 

0.08B£«a8e 

0.0808+888 

8.0088+000 

0.000EiB80 

8.8088+800 

OTTS/PHOTON 

RftNGE=245  11 

DEVIATION 

GROUP  LIBRARY  C 

MEASURED  C-10B0 

FROM  MEAN 

133e 

e.eeeE>0ee 

0.0008*000 

0.0008*000 

0.0808*000 

1000 

3.84EE-002 

2.7618-002 

'i.757E-a03 

0.0008*000 

2.2358*000 

800 

2.4')8E-002 

2.9428-002 

4.4378-003 

0.0008*000 

-.9588*000 

700 

1.527E-002 

I. 6678-002 

A.092E-003 

0.0008*000 

3.9048-001 

600 

1.655E-002 

1.5958-002 

3.3B8E-003 

0.0008*000 

1,7758-001 

512 

5.258E-00A 

0.0008*000 

1.03AE-003 

0.0008*000 

5,0978-001 

510 

1.AA7E-802 

1.7898-002 

4.765E-003 

0.0008*000 

-.7158*000 

450 

1.3078-002 

1.6518-082 

5.A0AE-003 

400 

3.312E-002 

0.0008*000 

-.6378*000 

3.3828-002 

9.305E-003 

0.0008*000 

-.7118-001 

300 

4.237E-002 

4.7308-002 

1.702E-002 

0.0008*000 

-.2898*000 

200 

6.5768-032 

1.0058-001 

1.4838-002 

1  <^0 

0.0008*000 

-.2348*001 

t  Ub 

!.:i 78-001 

1.2128-001 

2.0418-002 

0.0008*000 

2.4248-002 

100 

9.5318-002 

3.5748-002 

2.0998-002 

0.0008*000 

4. 5588-001 

""  8.1318-002 

9.3098-802 

1.3368-002 

0.0008*000 

-.6118*000 

60 

9.1458-002 

1.0658-081 

2.5328-002 

0.0008*000 

-.5948+000 

"  4.7208-002 

2.7918-002 

2.2498-002 

0.0008*000 

0,5818-001 

3.0658-003 

5.7138-003 

2.4078-003 

0.0008*000 

-.1108*001 

20 

7.6188-005 

0,0008*000 

2.7708-001 

2.7518-004 

0.0008*000 

10 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 
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COUNTS/PHOTON 

RAN6Es333  H 

DEVIATION 

GROUP  LIBRARY  C 
1330 

0.0eBE^000 

e.eeeE^eea 

HEASURED  0-1000 

FROM  MEAN 

0. 0000^000 

0.8000^000 

0.0000*000 

1000 

1.979E-002 

2.609E-003 

800 

1.271E-002 

2.723E-003 

1.561E-0B2 

0.030EV000 

1.6020*000 

9.3UE-003 

0.0000*000 

1.2490*000 

iw 

1.1900-002 

3.181E-003 

600 

l.lBlE-002 

2.934E-003 

1.270E-002 

0.0000*000 

-.2510+000 

9.182E-0B3 

0.0080*000 

8.9400-001 

5'2 

2.347E-004 

5.322E-004 

2.0490-004 

0.0000*000 

5.6000-002 

510 

4.407E-003 

3.681E-003 

450 

8.504E-003 

4.224E-003 

***  2.312E-002 
7.014E-003 

6.3500-003 

0.0000*000 

8,3050-001 

6.1160-003 

0.0000*000 

5.6510-001 

1.8820-002 

0.0000*000 

6,1410-001 

300 

3.213E-0I2 

3.920E-B03 

4.6310-002 

0.0000*000 

-.1590*001 

200 

3.632E-B02 

1.036E-002 

3.9270-002 

0.0000*000 

-.2840+000 

150 

7.143E-B02 

1.523E-002 

100 

i.537E-002 

2.301E-002 

5.410E-002 

1.952E-002 

9.2000-002 

0.0000*000 

-.6930*000 

5.2170-002 

0.0000*000 

5.5460-001 

3.6580-002 

0.0000*000 

8.9710-001 

60 

6.494E-002 

2.269E-002 

45 

3.618E-002 

1.0BSE-002 

7.6720-002 

0.0000*000 

-.5190*000 

4.0480-002 

0.0000*000 

-.3960*000 

30 

2.814E-383 

1.616E-003 

20 

S.960E-006 

4.184E-005 

2.5500-003 

0.0000*000 

-.3320*000 

0.0000*000 

0.0000*000 

2.1420-001 

0.000E*080 

0.000E«000 

0.0000*000 

0.0000*000 

0.0000*000 

Ill 


COUMTS/PHOTON 
BROUP  LIBRftRY  C 
1338 

a.eesE^aee 

0.e88E*e8i 


1.195E-002 

1.352E-a0I 

9.126E-883 

1.874E-003 

7,472E-003 

2.886E-003 

8.631E-083 

3.149E-003 

3.052E-004 

4.752E-004 

7.473E-003 

2.430E-003 

7.598E-003 

2.481E-003 

1.4S6E-e02 

4.576E-003 

2.135E-802 

4.750E-003 

2.109E-002 

4.6488-003 

4.690E-002 

1.020E-002 

3,8168-002 

1.012E-002 

3.300E-002 

9.596E-003 

4.394E-002 

l.:45E-002 

'  2.482E-002 
8.229E-B03 

I 

1.676E-003 

1.602E-003 

4.389E-006 

1,960E-00S 

I 

8. 0008+000 
0.000E+000 


RAN6E«39e  H 

flEASURED  8-1000 

0.0008+000 

0.0008+000 

1.2428-002 

0.0008+000 

1.3008-002 

0.0008+000 

9.3528-003 

0.0008+000 

7.4538-003 

0.0008+000 

2.1088-005 

0.0008+000 

9.7218-003 

0.0008+000 

9.0418-003 

0.0008+000 

1,3198-002 

0.0008+000 

1.7528-002 

8.0008+800 

1.3648-002 

0.0008+000 

5.6718-002 

0.0008+000 

4.0688-002 

0.0008+000 


4.402E-002 


2.3908-002 

0.0008+000 

3.8598-004 

0.000E+000 

0.000E+000 

0.0008+000 

0.0008+000 

0,000E+000 


DEVIATION 
FROn  NEAN 


0.0008+000 


-.3448+000 


.2068+001 


-.6518+000 


3.7178-001 


5.9798-001 


-.9248+000 


-.5818+000 


3.6448-001 


8.0758-001 


1.6038+000 


-.9628+000 


-.2498+000 


4,4938-002 


-.631E-002 


1.1178-001 


8.0568-001 


2.2398-001 


0.8008+000 


CQUNTS/FHQTCN 

RAN6E=2B.  H 

DEVIATION 

GROUP  LIBRARV  A 
133e 

a.0eeE*eeR 

I1EASURED  A-iB0 

FROM  MEAN 

B.BBBE^BBB 

0.00BE+0e0 

e.eeeE^eea 

B.0BBE9080 

leee 

5.':40E-002 

4.66BE-B82 

1.3i5E-003 

0.8BBE4-BBB 

9.383E+888 

900 

3.0A3E-002 

1.541E-0B2 

1.140£-0a3 

B.BBBE^BBB 

i.318E+88i 

700 

2.547E-001 

2.528E-0B1 

3.:T.3E-003 

B.BBBE^BBB 

1.297E+001 

600 

i.l6lE-002 

0.000E9000 

1.0926-003 

0.B00E*000 

1.0640+801 

5  i  2 

l,811E-00.i 

1.539E-001 

1.479E-002 

B.BBBE^BBB 

1.8350+008 

510 

1 . 96oE'002 

0.008E90B0 

1 . 138£*003 

0.0B0Et000 

1.64BE+8BI 

450 

1.472E-0B2 

5.716E-B03 

1.323E+BB1 

l,437E-803 

0.800E9008 

400 

7.0690-002 

1.087E-001 

-.151E+802 

2.504E-003 

0.B8BE9B80 

300 

6.154E-002 

7.729E-083 

2.760E-003 

B.BBBE^BBB 

1.950E+081 

200 

1.183E-001 

1.491E-B01 

-.957E+001 

3.221E-003 

0.BBBE9B8B 

1  JU 

9.528E-O02 

1.094E-0B1 

-.47BE+8B1 

2.9B9E-003 

B.BBBE^Bee 

100 

1.492E-002 

2.235E-002 

-.466E+001 

1.593E-003 

e.BBBE^BBe 

‘  2.083E-002 

i.309E-0B2 

7.767E+8B0 

7.971E-004 

B.BBBEiBea 

60 

2.531E-002 

3.929E-002 

-.1B4E+882 

1.339E-003 

0.8BBE9000 

JTj 

1.524E-002 

3.5070-803 

1.164E-003 

30 

3. £9 /E -004 

0.eeeE^003 

1.0070+881 

0.0080^000 

i .652E-004 

0.8080^000 

2.359E+000 

20 

0.0000^000 

0.008Ei000 

0.000E+000 

O.000E^000 

0.0080+880 

10 

0.200E+000 

0.8000+000 

0.000E*000 

0,8000+000 

0.0000+000 
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COUNTS/PHOTOM 

RANGEMl  K 

DEVIATION 

GROUP 

1330 

LIBRARY  A 

MEASURED  A-10B 

FROM  MEAN 

I.BBBE4B0B 

8.B0IE4B0B 

8.000E*000 

1  ama 

0,B0eE*B00 

8.Be0E*000 

II0I0V 

5.024E-002 

1.441E-003 

^'0?§00E!iB8 

-.706E*0e0 

B00 

2.6B4E-002 

1.473E-B02 

1.1678*001 

1.B37E-003 

0.000E*B00 

700 

2.449E-00i 

2.475E-001 

-.1088*001 

2.42SE-003 

0.800E*000 

600 

2.253E-002 

4. 4488-002 

-.1508*002 

1.454E-003 

0.000E*000 

512 

1.437E-001 

1.149E-001 

1.5S3E*081 

1.85IE-003 

0.000E*000 

510 

3.117E-002 

5.312E-002 

-.1498*002 

450 

1.474E-003 

0.000E*080 

3,306E-002 

5.969E-002 

-.1828*002 

400 

1.458E-003 

B.000E*000 

a.656E*002 

7,9048-802 

2.7078*000 

300 

2.778E-003 

0.0088*000 

9,43e£-002 

4.971E-002 

1.5098*001 

2.961E-003 

0.0008*000 

200 

l.20SE-B0i 

1.407E-001 

-.621E*001 

3.251E-003 

0,000E*000 

150 

1.200E*801 

1.233E-001 

-.8388*000 

1  aa 

3,963E-003 

0.0088*800 

4.043E-002 

1.1288-002 

9.9848*000 

75 

2.919E-803 

0.0008*000 

3.013E-002 

5.1408-003 

1.2508*001 

60 

1.998E-003 

0.0008*000 

2.934E-002 

0.0088*000 

1.945E-003 

0.0008*000 

1.5098*001 

1.580E-002 

0.0008*000 

1.4248*001 

1.109E-003 

0.0008*000 

30 

7.779E-B04 

0.0008*000 

5, 1928*000 

1.498E-004 

0.0008*080 

20 

1.608E-005 

0.0008*800 

1.3118*000 

1.227E-005 

0.0008*000 

10 

0.000E*000 

0.0008*000 

0.0008*800 

0.0008*000 

0.0008*000 

COUNTS/PHOTON 

RAN6E>83  H 

DEVIATION 

6R0UP  LIBRARY  A 

MEASURED  A- 100 

FROM  HEAR 

1330 

8.e0eE^000 

0.ee0E^000 

0.0005*000 

0.0005*000 

0.0005*000 

1000 

3.558E'002 

9.961E-004 

2.6655-002 

0.0005*000 

8.9625*000 

800 

2.122E-e02 

8.7B9E-004 

5.2355-002 

0.0005*000 

-.3545*002 

7Ra 

1.677E-001 

2.040E-003 

1.0225-001 

0.0005*000 

3.2125*001 

600 

3.517E-002 

l.a4SE-003 

4.3925-002 

0.0005*000 

-.4735*001 

512 

B.760E-002 

1.424E-003 

1.5655-001 

8.0005*000 

-.4835*002 

510 

3.923E-002 

1.714E-003 

4.3595-002 

0.0005*000 

-.2545*001 

450 

4.0005-002 

1.897E-003 

6.7305-003 

0.0005*000 

1.7525*001 

400 

1.021E-001 

3.537E-003 

8.6085-002 

0.0005*000 

4.5315*000 

300 

1.260E-001 

3.851E-003 

9.2165-002 

0.0005*000 

8.7905*000 

200 

1.198E-001 

3.652E-003 

2.1505-001 

0.0005*000 

-.2605*002 

150 

1.6445-001 

4.377E-003 

1.7205-001 

0.0005*000 

-.1755*001 

100 

7.309E-002 

3.828E-B03 

75 

4.653E-002 

2.598E-003 

0.0005*000 

0.0005*000 

1.9095*001 

4.1655-002 

0.0005*000 

1.8305*000 

60 

4.580E-002 

2.780E-003 

0.0005*000 

0.0005*000 

1.6475*001 

45 

2.174E-002 

2.126E-003 

30 

1.177E-003 

2.379E-004 

8.8945-003 

0.0005*000 

6.0415*000 

8. 8345-004 
0.0005*000 

1.2345*000 

20 

0.0005*000 

0.0005*000 

0.0005*000 

0.0005*000 

0.0005*000 

10 

0.0005*000 

0.0005*000 

0.0005*000 

0.0005*000 

0.0005*000 

I 
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I* 


GROUP 

1330 

LIBRARY  A 

HEASURED  A-i( 

0.0B0EtBB8 

0.0008*000 

laaa 

0.0008^000 

0.0008*000 

2.633E-B02 

2.414E-002 

8.9878-004 

a.000E*000 

aaa 

1.S73E-B02 

8.935E-0B3 

700 

6.613E-004 

0.0008*000 

1.1B1E-0B1 

1.B14E-B01 

600 

1.869E-003 

0.0008*000 

3.2BSE-002 

3.2188-002 

512 

1.219E-0a3 

0.0008*000 

5.572E-002 

5.375E-002 

510 

8.090E-004 

0.0008*000 

3.796E-002 

4.4808-082 

450 

1.477E-003 

0.0008*000 

3.0958-002 

3.674E-002 

1.415E-003 

0.0008*000 

400 

9.565E-002 

6.9608-002 

300 

2.786E-003 

0.0008*000 

1.350E-001 

1.3018-001 

4.057E-003 

0,0008*000 

200 

150 

l.i63E-081 

3.652E-003 

1.2888-001 

0.0008*000 

3.9628-002 

1.704E-001 

100 

6.046E-003 

0.0008*000 

1.082E-001 

1. 3218-001 

75 

S.SI3E-003 

0.0008*080 

7.551E-002 

1.1018-081 

60 

3,734E-003 

0.0008*008 

6.651E-0B2 

6.5238-002 

45 

4.221E-003 

0.0008*000 

3.209E-002 

1.8548-002 

30 

2.477E-003 

0.0088*000 

2.036E-003 

2.6388-805 

20 

3.965E-004 

0.0008*000 

3.814E-006 

0.0008*080 

10 

2.718E-006 

0.0008*000 

0. 0008*000 

0.0008*080 

0.0008*000 

0,0008*000 

DEVIATION 
FRON  NEAN 

a.aeflE^Bea 

2.443E^aae 

i.a28Eiaai 

8.953E+aaa 

5.444E-Bai 

2.435EiaaB 

•.4&2E+Bai 

4.e91E^B8B 

9.352E+8BB 

1.189E*Bee 

-.341E+eBl 

2.i63Etaei 

-.433E+BB1 

-.925E+BBI 

3.a31E-Bei 

5.473E+eBB 

s.a68EteBe 

i.4e3Eteee 

e.eeeE^BBe 
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COUNTS/PHOTON  RAN6EM6&  H 


DEVIATION 


SROUP 

LIBRARY  A 

NEASUREO  A-t00 

FROM  KEAN 

1330 

8.B0BE*e0B 

0.0000*000 

0.0000*000 

0.0000*090 

0.9000*009 

1000 

1.914E-002 

a.086E-9B4 

1.8780-002 

0.0090*900 

2.9100*000 

800 

1.1270-002 

7.047E-004 

7.3460-003 

0.0000*009 

5.5670*090 

700 


8.1230-002 

7.3570-002 

1.2380-003 

0.9000*990 

6.1870*999 

600 

2.8800-002 

3.5460-992 

1.0530-003 

0.0000*009 

-.6320*001 

512 

3.5550-002 

6.9160-002 

6.0270-004 

0.0990*009 

-.5570*002 

510 

3.5140-092 

1.7740-002 

1.0020-003 

0.0000*000 

1.7370*001 

4se 


490 


380 


3.432E-O02 

1.30SE-093 

8.ai6E*002 

2.741E-003 


1.781E-0e2 

9.oeeE^0ao 

1.230E-00i 

0.000E*0aa 


i.327Eiaei 

-.127E+002 


3.8980-003 

0.0090*000 

1.7140*008 

200 

9.9550-002 

5.7160-002 

3.9430-003 

0.0900*900 

1.0750*001 

150 

1.6350-001 

1.9290-001 

4.7870-003 

0.0900*000 

-.5950*001 

100 

1.1720-001 

3.1480-002 

5.9760-003 

0.0800*000 

1.4350*001 

75 

7.6790-902 

2.8220-902 

3.6690-003 

0.0000*000 

1.3230*001 

60 

8.7650-002 

5.4480-902 

45 

5.3650-003 

0.0900*000 

6.1830*000 

4.2600-002 

1.1560-002 

3.5290-003 

0.0900*000 

8.7970*000 

30 

2.9120-003 

8.7990-003 

3.7450-004 

0.0000*000 

-.1570*002 

20 

1,7390-006 

8.8000*000 

9. 9310-007 

0.0000*000 

1.7510*000 

10 

0.0090*000 

0.0000*000 

0.0090*000 

0.0000*000 

0.0000*000 
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viii:'Vt=;F«OTrK 

DEVIATION 

GROUP  LlBRftRV  A 

HEASUREO  A- 100 

FROM  HEAN 

1  ^  50 

i.eeeEmi 

o.aoeE^fles 

0. 0000*000 

0.000E*000 

0.0008*000 

leeo 

3.929E-a84 

aea 

6.137E-0a3 

3.69iE-a04 

708 

4.0e7E-002 

'i,661E-004 

1.09BE-002 

0,000E*000 

-.6098*000 

9,921E-003 

0.000E*000 

-.1028*002 

3.466E-002 

0.000E*000 

5.5958*000 

600 

2.041E-002 

8,336£-004 

3.192E-002 

0. 0008^000 

-.1388*002 

512 

1.438E-002 

3.573E-004 

510 

2.:52E-002 

3.902E-004 

7.156E-003 

0.000E*000 

2.0208*001 

1.08iE-002 

0.000E*000 

1.4278*001 

450 

2.371£-aB2 

a.396£-004 

1.774E-002 

0.0008 *080 

7.1068*000 

400 

6.089E-002 

2.239E-a03 

300 

9.273E-a02 

3.549E-003 

:,382E-002 

0.0008*000 

1.6568*001 

4.366E-002 

0.0008*000 

1.3828*001 

■>00 

.4B2E-002 

l,851E-003 

7.793E-002 

0.0008*000 

-.7028*001 

150 

1.292E-001 

4,16aE-003 

7.613E-082 

0.0008*000 

1.2748*001 

100 

9.034E-002 

3.568E-003 

■yc 

1.0028-001 

0.0008*000 

-.2628*001 

■'  7.496E-a02 
7.514E-003 

2.235E-002 

0.0008*000 

7.0018*000 

8.839E-002 

5.567E-003 

:. 5 198-003 
0.0008*000 

1.5438*001 

Jf 

■’  4.720E-002 
i.o85E-003 

*  11 

3.7478-002 

0.0008*000 

3,6238*000 

4.i;4E-0a3. 

9.346E-004 

5.7178-005 

0.0008*000 

4,4058*000 

"0 

l.:34E-006 

7.773E-00/ 

0.0008*000 

0.0008*000 

1.5378*000 

1  0 

0.080E^000 

0.000E»000 

0.0008*000 

0.0008*000 

C.000E*000 

COUNTS/PHQTON 

RAN6E:333  N 

DEVIATION 

GROUP  LIBRARY  A 

MEASURED  A-100 

FROM  MEAN 

1338 

e.esaE^eefl 

e.eflflEteee 

0.00084000 

0.00884000 

0.00084008 

S.488E'883 

2.168E-884 

OAil 

7.3238-003 

0.00084000 

-.84684001 

OvIq 

3.5038-883 

2.155E-004 

'inn 

4.6368-003 

0.00084080 

-.52584001 

1.985E-082 

4.154E-084 

Ltlll 

2.5398-002 

0.00084000 

-.13384002 

Oov 

1.132E-002 

4.6248-884 

512 

5.4888-883 

2.4918-004 

510 

1.4858-082 

5.9778-004 

1.3468-002 

6.7698-004 

3.9598-002 

9.6718-804 

5.3238-003  , 
0.00084000 

1.29684001 

6.3028-003 

0.00084000 

-.33084001 

2.4428-003 

0.00084080 

2.07684081 

6.0688-083 

0.00884008 

1.09184001 

1.4458-082 

0.00084000 

2.60084001 

380 

6.3888-002 

5.6778-003 

7.8578-002 

0.00084000 

-.25884001 

200 

4.2488-002 

1.8068-003 

4  en 

3.1378-002 

0.00084000 

6.15384080 

1  wv 

8.3368-002 

3.1028-003 

4  na 

9.1688-002 

0.00084080 

-.26884001 

1  vV 

6.4558-082 

2.4718-003 

•jg 

7.8998-002 

0.80084000 

-.58484001 

6.  1478-002 
4.2018-003 

5.1748-002 

0.00084000 

2.31584000 

60 

6.4088-002 

3.7568-003 

1.0908-001 

0.00084000 

-.11984002 

3.1868-002 

2.1778-003 

1.0118-002 

0.00084080 

9.99384000 

2.3768-003 

4.5048-004 

7.8638-004 

0.00884000 

3.53084000 

id 

5.4738-006 

3.3068-006 

4  n 

0.00084000 

0.00084000 

1.65684000 

lu 

0. 0008^000 
0.000Et000 

0.00084000 

0.00084000 

0.00084008 
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COUNTS/PHOTON  RftN6E=20.  N 


DEVIATION 


GROUP 

'  LIBRARY  B 

K8A3URED  6-100 

FROH  H8AN 

1330 

0.000E>000 

0.000E>000 

0.0008+000 

0.0008+000 

0.0008+000 

1000 

0.000E*000 

0.000E^000 

0.0008+000 

0.0008+000 

0.0008+000 

800 

0.000E^000 

0.000EV000 

0.0008+000 

0.0008+000 

0.0008+000 

700 

0.000E^000 

0.000E^000 

0.0008+000 

0.0008+000 

0.0008+000 

800 

512 

0.000E^000 

0.000E^000 

0.0008+000 

0.0008+000 

0.0008+000 

0.000E^000 

0.0008+000 

0.000Et000 

0.0008+000 

0.0008+000 

510 

0.000E^000 

0.0008^000 

0.0008+000 

0.0008+000 

0.0008+000 

A50 

7.940E-0B2 

1.832E-003 

7.5788-002 

0.0008+000 

1.9818+000 

400 

i.463E-002 

1.930E-003 

4.3348-002 

0.0008+000 

2.8578+001 

300 

2.445E-002 

2.052E-003 

0.0008+000 

0.0008+000 

1.1928+001 

200 

8.608E-'001 

3.061E-003 

8.2618-001 

0.0008+000 

1.1288+001 

150 

1.759E-001 

7.826E-003 

1.5938-001 

0.0008+000 

2.1758+000 

100 

T5 

3.i;4E-B02 

2.743E-003 

2.5858-002 

0.0008+000 

2.1478+000 

1. 8798-003 

0.0008+000 

1.2398-003 

0.0008+000 

3.9378+000 

;. 6808-003 

0.0008+000 

45 

3.0978-004 

0.0008+000 

2.0518+000 

0.e008<'000 

0.0008+000 

’  I 

0.0008^000 

0.0008+000 

0.0008+000 

0.e00E»000 

0.0008+000 

20 

0.0008^000 

0.0008+000 

0.0008+000 

0.0028>000 

0.0008+000 

1  i"? 

0.0008^000 

0.0008+000 

0.0008+000 

lu 

0.0008t000 

e. 0008+000 

0.0008+000 

0.0008+000 

0.0008+000 
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i  !• 

i 

I 


COUNTS/PHOTON 
GROUP  LIBRARY  B 
1330 

B.eBBE^BBB 

0.88eE*8e8 


8.e8eE>0B8 

e.8eeE^a8a 


e.aaa^aaa 

8.880E*a8a 

a.aeeE^eae 

B.eeeE^eee 

a.aeaE^aae 

8.88aE^888 

а. eaaE^aae 

e.aaaE^aaa 

e.aaaE^aea 

6.4fa8E-a02 

1.477E-a03 

'  7.946E-aa2 
2.8i&E-ee3 

'  4.089E-002 
2.4<i7E-803 

'  5.425E-001 
4.698E'ee3 

'  2.726E*0ai 

б. 847E-ee3 

1 

7.534E-0a2 

4.332E-003 

'  3.223E-a02 
4.t23E-ae3 

°  i.e45E-0a2 
2.a9aE-8e3 

'  2.8ME-083 
7.l84E-aa4 

^  7.984E-085 
6.8988-885 

*  3.091E-006 
3.09lE-ea6 

^  8.8888^888 


RAN6EM1  H 

HEASURED  8-188 

DEVIATION 
FRON  HEAN 

0. 0088*808 

B.e08E*080 

8.0808*880 

8.0008*808 

0.8a0E*aaa 

0.0008*800 

8.8001*000 

8.8a0E*a0a 

0.0088*008 

8.8008*000 

8.8808*800 

8.8888*000 

8.8888*888 

8.8808*888 

8.0808*080 

8.8888*000 

0.0088*008 

8.8008*888 

8.8808*800 

8,0808*800 

0.8888*800 

1.3438-002 

0,0808*000 

3.4658*801 

5.2968-002 

8.8808*880 

1.3158*801 

1. 8488-083 
8.0008*808 

1.5588*081 

5.2678-001 

8.0088*880 

3.3608*800 

5.3258-001 

0.8808*808 

-.4298*002 

1.0578-001 

8.0008*800 

-,7008*801 

8.0008*000 

8.0008*000 

7.8158*008 

8.0088*088 

0.0808*800 

4.9978*000 

8.8008*080 

8.8008*800 

2.8688*880 

8.8808*880 

8.0888*080 

1,3118*080 

8.8808*000 

8.0888*888 

1.0008*008 

8.8808*088 

8.8008*808 

8.0808*008 

a 


121 


I 


CO'JNTS/PHOTON 

RANGE:83  H 

DEVIATION 

GROUP  LIBfiflfiV  fl 

MEASURED  B-100 

FROM  MEAN 

1330 

e.BBeE^eee 

0.000E>000 

e.000E^000 

0.000E*000 

0.0008*000 

1000 

0.000E^000 

0.000E*000 

ann 

0.000E>000 

0.000E*000 

0.0008*000 

aou 

1.651E-005 

1.651E-00S 

700 

0.000£*000 

0.000E*000 

0.000E40B0 

0,0e0E*000 

1.0008*000 

0.O00E«000 

0,000E^000 

0.0008*000 

j)00 

0.0e0E^000 

0.000E*000 

0.e00E40e0 

0.00eEt000 

0.0008*000 

513 

0.000E^000 

0.000E>0e0 

510 

0.000EV000 

0.000E*000 

0.000E>000 

0.000E4000 

0.0008*000 

0.000E>000 

0.008E<-000 

0.0008*000 

AS0 

4.199E-002 

1.U0E'003 

400 

6.441E-002 

1.438E-003 

300 

4,505E-002 

2,153E-003 

3,922E-002 

0,000E^000 

2.497E*000 

7.507E-002 

0.0e0E*000 

-.7418*001 

3.570E-002 

0,000E*0e0 

4.3428*000 

200 

3,559E-001 

3.248E-003 

3.068E-001 

0.000E*000 

1. 5118*001 

150 

3.408£-e0l 

4.401E-003 

1.847E-001 

0.000E*000 

3.5468*001 

100 

1.62<iE-001 

a.047E-003 

9.101E-002 

0.000E*000 

1.1838*001 

3..=!-7E-002 

-,43l£-003 

3.339E-002 

0.000E^000 

3.0718*000 

:.a:.0E-002 

3.755E-003 

0.000E>000 

0. 0008*000 

1,4998*001 

l,7l8£-002 

2.140E-003 

0.0008*000 

0.0008*000 

8.0298*000 

t'./'52E-004 

2.596E-004 

0.0008*000 

0.0008*000 

3.3878*000 

0.000E>000 

0.000E+008 

0.0008*000 

0.0008*000 

0.0008*008 

0.000E^00e 

0.000Et000 

2. 0008*000 
0.0008*000 

0.0008*000 

I 


i  %• 


COUNTS/PHOTON 
GROUP  LIBRARY  B 
1330 

0.eBIE+000 

0.eB8Et000 

1000 

0.e00E^000 

0.00BE*0B0 

000 

0.008E^0B0 
0.000EV000 


700 


600 


512 


510 


450 


400 


300 


200 


150 


100 


60 


30 


20 


10 


..000E>000 

0.000E*B00 

0.000E>000 

0.000E*0B0 

0.000E>0B0 
.000E>000 

0.000E^008 

0.B00E>000 

2.70<JE-002 

6.676E-004 

4.565E-002 

1.472E-003 

4.784E-0B2 

2.829E-003 

2.282E-801 

2.9B0E-003 

3.164E-001 

3.215E-003 

l.862E'eBl 

8.803E-003 

1.167E-001 

5.147E-B03 

1.144E-001 

8.127E-003 

3.946E-002 

2.787E-003 

3.158E-003 
5.054E-004 

1.382E-005 

9.368E-006 

0.000E^000 

0.000E«’880 


RAN8EM24.S  H 
MEASURED  B-108 

0.000E4BB0 

0.000E^800 

0.008Ei00e 

0,000E^008 

0.008E^008 

0.000E>000 

0.000E+000 

0.008Et000 

0.0B0E>000 

0.000E>000 

0.08BE4BB0 

0.00eE^000 

0.000Ei0B0 

0.00BE^0B0 


5.498E-0B2 

0.008E^000 


0.001 


>000 


3.025E*001 

0.000E^000 

2.746E*001 

0.000^000 

1.824E-001 

0,000E*000 


5.54 


.  J02 

I0E*000 


5.009E-002 

0.B00E^000 

0.000Et000 

0.000E^800 

5.358E-004 

0.000E*000 

0.000Et000 

0.000E^080 


,0008*000 


DEVIATION 
FROM  MEAN 


0,000E4000 

0.000E4000 

0.000E*000 

0.000E*000 

0,000E*000 

0,000E*000 

0,000E*000 

-.25IE+082 

•.633E4001 

-.7l8Et001 

-.255E*002 

5.B82E*000 

4.282E-001 

1.189E*001 

7.909E+000 

1,416E*001 

4.479E+000 

1.476E*000 

0.000E*000 
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COUNTS/PHOTQN 

RANGE-166  « 

DEVIATION 

GROUP  LIBRARY  B 
1330 

HEA3URE0  B-10B 

FROH  MEAN 

B.BBBE+Bee 

0.0003*000 

0.0003*000 

0.0e0E>ee0 

0.0003*000 

1 ODV 

B.eeeE^eoe 

0.0003*000 

0.0003*000 

0.0eeE4000 

0.0003*000 

300 

0.030E*00B 

0.0003*000 

0.0003*000 

0.000Et000 

0.0003*000 

700 

0.000Et000 

0.0003*000 

e. 0003^000 

■  HA 

0.0003*000 

0.0003*000 

QOld 

■.02YE-005 

8.0003*000 

1.0003*000 

7.02')E-005 

0.0003*000 

0.030E*000 

2.0003*000 

0.0003*000 

0,0003^000 

5i0 

0.O003>000 

0.0003*000 

2.2003*000 

0. 0003^000 

0.8003*000 

0.0083*000 

430 

1.6213-002 

2.8863-002 

-.2743*002 

4.6063-004 

0.0003*000 

400 

3.568E-002 

3.1753-003 

2.5133*001 

1.2933-003 

0.0003*000 

300 

4.356E-002 

4,0363-002 

1.6433*000 

1.9443-003 

0.0003*000 

1.4923-001 

2.5423-003 

150 

2.4273-001 

1.4843-001 

0.0003*000 

3.4583*000 

2.5673-001 

-.2833*001 

4.9613-003 

0.0003*000 

100 

1.7573-001 

3.0963-001 

-.2173*002 

6.1443-003 

0.0003*000 

1.276E-001 

2.2213-001 

-.1373*002 

5.120E-003 

0.0003*000 

cl* 

'  1.2:93-801 

5,6753-002 

1.0013*001 

5.9113-003 

0.0003*000 

■  5.9493-002 

1 .  .'623-002 

6.4193*000 

3.4073-003 

1  A 

0.0003*000 

316E*233 

6.6123-004 

3.4413*000 

4.33aE-«04 

0.0003*000 

-  2 

;. 7023-006 

2.0003*000 

1.2203*000 

1,2143-006 

0.2003*000 

'  0 

0.300Et000 

2.0003+000 

0.0003*000 

0. 0003*000 

0.0003*000 
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COUNTS/PHQTON 

RANGE=245  H 

DEVIATION 

GROUP  LIBRARY  B 

MEASURED  B-100 

FROH  HEAN 

1330 

1.708E-005 

1.700E-005 

1000 

0.000E*800 

0.e00E^000 

0.000E^000 

0.000E*000 

1,000E^000 

0.000E*000 

0.000E>000 

0.000E*000 

Oou 

0.030E*000 

0.000E*000 

0.e00E^000 

0.000E^000 

0.0008*000 

700 

0.00eE*000 

0.000EV000 

L  fin 

O.000E4000 

0.000E^000 

0.00BE*000 

ObD 

0.000E*000 

0.000E^000 

0.000E*000 

0.000Ef000 

0.0008*000 

c  1 0 

0.000E^000 

0.000E*000 

0.008E^000 

0.000E>080 

0,000E*000 

510 

0.000E>000 

0,000Et000 

0.000E*080 

0.000E^000 

0.000E*000 

450 

6.a80E-0e3 

3.294E-004 

6.317E-003 

0.000Et000 

1.710E*000 

400 

1.753E-002 

S.926E-004 

9.097E-B03 

0.000E^000 

1.4248*001 

’00 

'  2.650E-002 

1.191E-003 

200 

6.94e£-002 

2.023E-003 

2.176E-002 

0.0eaE>000 

3.900E*000 

5.6a8E-002 

0.008E*000 

6,626E*000 

1.379E-001 

3.260E-003 

100 

l.lllE-001 

3.399E-003 

75 

9.930E-002 

3.835E-003 

60 

1.124E-001 

4.555E-003 

6.524E-002 

0.000E*000 

2.229E*00l 

1J33E-001 

0.000Ef000 

-.182E*002 

8.860E-002 

0. 0008^000 

2.7918*000 

3.062E-001 

0.00BE>000 

-.4258*002 

4^ 

■'  5.875E-002 
2.441E-003 

7.503E-002 

o.eoBE^-eBe 

-.616E*001 

'  21 

3.274E-003 

4.004E-004 

20 

1 . 766E-005 
1.177E-005 

3.304E-004 

0.0008^000 

7.3528*000 

0.000E*00e 

0.000E>000 

1.501E*000 

1  K 

0.000E^300 

0.000E>0e0 

0.000E^000 

B.000E»000 

0.0008*000 

f 
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COUNTS/PHOION 

RANGE=333  « 

DEVIATION 

GROUP  LIBRflRy  B 

MEASURED  B-100 

FROM  MEAN 

1330 

0.000E^00e 

0.0008*000 

0.000E^000 

0.0008*000 

0.0008*000 

1000 

0.000Et000 

0.0008*000 

0.000E*000 

0.0008*000 

0.0008*000 

800 

0.000Ei000 

0.0008*000 

0.000Et000 

0.0008*000 

0.0008*000 

700 

0.000E^000 

0.0008*000 

0.000E^000 

0.0008*000 

0.0008*000 

&00 

0.0e0Ei000 

3.0008*000 

0.000E>000 

0.0008*000 

0.0008*000 

c  ^2 

0.000E^000 

0.0008*000 

0.000E>000 

0.0008*000 

0.0008*000 

510 

0.000E>e00 

0.0008*000 

0.000E^000 

0.0008*000 

0.0008*000 

450 

':.635E-003 

0.0008*000 

2.080E-004 

0.0008*000 

1.2678*001 

400 

7.854E-003 

0.0008*000 

3.485E-004 

0.0008*000 

2.2658*001 

300 

1.467E-002 

1.7288-002 

9.581E-004 

200 

2.836E-002 

0.0008*000 

-.2518*001 

2.7658-002 

1.262E-003 

1 

0.0008*000 

5.6028-001 

i.886E-002 

4.9718-002 

2.645E-003 

0.0008*000 

7.2428*000 

100 

5.934E-082 

9.8758-002 

3.433E-003 

0.0008*000 

-.1148*002 

75 

4.870E-002 

6.5088-002 

2.237E-003 

0.0008*000 

-.7328*001 

60 

6.9968-002 

9.7618-002 

3.705E-003 

0.0008*000 

-.7468*001 

‘  4.081E-002 

6.1758-002 

4.777E-003 

0.0008*000 

-.4388*001 

30 

3.0078-003 

5.0168-005 

4.638E-004 

"0 

6.250E-006 

0.0008*000 

6.3748*000 

0.0008*000 

0.0008*000 

4.078E-006 

10 

0.0008^000 

1.5338*000 

3.0008*000 

0. 0008^000 

0.0008*000 

0.0008*000 
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■ 

■ 

P. 


i  If 

i 


COUNTS/PHOTON 
GROUP  LIBRARY  C 
1338 

0.eeeE>e8t 

8.8880^088 

1008 

2.285E-001 

2.&4&E*883 


1.1210-801 

1.779E-003 

3.0110-082 
1.8160-083 

8.7950-003 
8.6310-084 

4.4170-004 

1.9020-004 

6.8380-883 

7.7950-084 

5.9110-083 

6.9140-004 

1.2960-082 

1.4740-003 

2.1100-002 

1.6400-003 

4.1680-001 

3.3180-003 

1.5360-001 

3.9630-003 

4.6890-002 
2.2410-003 


700 


600 


512 


510 


450 


400 


300 


200 


150 


75 


60 


45 


3.6530-003 

1.3160-003 

2.3270-003 

9.5760-004 

3.6380-004 

2.9670-004 

0.0000*030 


20 


10 


0.0000*000 

0.0000*000 

0.0080*000 

0.8000*000 


RAN8E=20.  N 

HEASURED  C-108 

DEVIATION 
PROn  NOAH 

0.0000*000 

0.0000*000 

0.0000*000 

^'0!0f0E*000 

-.7600*001 

2.8080-001 

2.6850-082  . 
0.0080*000 

3.2110*000 

7.8460-003  . 
0.0000*000 

1,0990*000 

«-n0 

2.3230*000 

8.2790-003 

0.0000*800 

-.1840*001 

**0?000!*S80 

8.5500*000 

^'0?000E*008 

-.3880*002 

*'r0ia0 

1.2860*001 

4.0310-001 

0.0000*000 

4.1330*000 

1.0090-001  . 
0.0000*000 

1.3320*001 

0.000E+00B_ 

0.0000*000 

2.0920*001 

2.0100-002 

0.0000*000 

-.1470*002 

0.0000*000 

0.0000*000 

2.4300*000 

0.0000*000 

0.0000*000 

1.2260*080 

0.0000*000 

0.0000*000 

0.0000*000 

0.0000*000 

0.0000*000 

0.0000*000 

0.0000*000 

0.0000*000 

0.0000*000 
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CCUNTS/PHOTON 
GROUP  LIBRARY  C 
133e 

e.BeaEioee 

e.aeeE^eae 


GBB 


700 


300 


1.832E-081 

2.509E-003 

';.833£-002 
1.618E-003 

3,3;5E-002 

1.196E-003 

1.312E-082 

1.129E-083 


510 


450 


2.563E-004 

1.171E-004 

1,0358-082 
1.0288-003 

8.699E-003 

9.8058-004 


2.909E-002 

2.0618-003 


200 


150 


100 


75 


60 


45 


30 


20 


3.1588-002 

2.4448-003 

3.4538-001 

3.8558-003 

2.289E-001 

5.036E-B03 

7.8058-002 

4.9r4E-003 

2.980E-002 

2.034E-003 

i.308£-002 

1.344E-003 

3.4^48-003 

I. 290E-003 

a.03iE-Z05 
5.  ?6iE-005 

0,000£t000 

0.000EV000 

0.000E^000 

J. 020E+000 


RANGE=41  H 

heasured  c-100 

0.0e0E*000 

0.0008^000 

1. 4098-001 
0.000Et000 

8.6258-002 

0.0008*000 

2.5028-002 

0.0008*000 

2.1348-002 

0.0008*000 

0.0008*000 

0.0008*000 

1.4998-002 

0.0008*000 

2.3108-002 

0.0008*000 

1.9068-002 


0.0008*000 

0.0008*000 

3.9128-001 

0.0008*000 

2,0418-001 

0.0008*000 

7.1008-002 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 

0.0008*000 


DEVIATION 
FROM  NEAN 


0.0008*000 

1.6868*001 

7.4848*000 

7.3018*000 

-.2858*001 

2.1838*000 

-.4518*001 

-.1468*002 

4.8678*000 

1.2928*001 

-.1198*002 


3.2898*000 

1.4118*000 

1.4658*001 

7.0918*000 

2.6628*000 

1.3438*000 

0.0008*000 

0.0008*000 
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CQUNT5/PHQT0N 

RAN6E=83  K 

DEVIATION 

BROUP 

1330 

LIBRARY  C 

HEASURED  C-i00 

FROH  HEAN 

0.000Et0B0 

0.0000*000 

0.0000*000 

f  Ma 

0. 0000*000 

0.0000*000 

10DV 

1.371E'001 

1.0160-001 

1,7520*001 

2.026E-003 

0.0000*000 

000 

7./32E-002 

7.4770-002 

1.5920*000 

708 
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Appendix  D.  Principle  of  Monte  Carlo 


Introduction 

The  Monte  Carlo  Method  uses  random  sampling  to  construct  the  solu¬ 
tion  of  a  physical  or  mathematical  problem.  This  random  sampling  dis¬ 
tinguishes  this  technique  from  numerical  techniques.  A  complete  treat¬ 
ment  of  the  Monte  Carlo  method  will  not  be  attempted  in  this  review, 
only  the  basic  principles  will  be  illustrated. 

Basic  Principles 

The  Monte  Carlo  method  depends  on  the  use  of  "pseudorandom"  numbers, 
?,  which  are  uniformly  distributed  on  the  Interval  (0,1),  i.e.  0<^C<^1. 

The  interested  reader  is  referred  to  several  references  that  discuss 
random  number  generation  (Refs  12:27-31;  3). 

The  basic  principle  is  based  on  modeling  a  physical  phenomenon  with 
the  proper  density  function.  If  the  phenomenon  can  be  modeled  as  n  in¬ 
dependent,  mutually  exclusive  events  Ej^,  E2.-.,E^,  with  probabilities 
Pj^,  P2»*-'>Pjj»  respectively,  and  ^  p^  =  1,  then  for  a  discrete  case 

E^  will  be  determined  by  the  generation  of  a  random  number,  i;,  on  the 
interval  (0,1).  P^  is  the  probability  density  function  and  is  illustra¬ 
ted  in  Figure  1. 

To  find  E^ 

P  i-1  "  Pi  ■'  2  Pi-1  i  Pi  P2---+  Pi  “Pi 

where  P(x)  is  the  cumulative  probability  distribution  function  for  the 
discrete  case.  Figure  2  illustrates  the  distribution  function. 


Finding  for  a  continuous  case  may  be  approached  in  a  similar 
fashion.  If  p(x)dx  is  the  probability  of  x  lying  between  x  and  x+dx, 
with  a£x£b,  and 

b 

/p(x)dx  *  1 
a 

then 

C  "P  (x)  =  Jpd 
a 

determines  x  uniquely  as  a  function  of  if  i;  is  uniformly  distributed 
on  Oj<c;<l,  then  x  falls  with  frequency  p(x)dx  in  the  Interval  (x,  x+dx). 

To  Illustrate  this  first  method  let  us  sample  the  distance  to  col¬ 
lision  of  a  particle.  The  density  function  for  this  example  is  given  by 

p(l)  =  S^e  (1) 

where 

1  =  Distance  to  collision 

=  Total  macroscopic  cross-section  of  the  medium.  is 
interpreted  as  the  probability  per  unit  length  of  a 
collision. 

The  probability  of  a  first  collision  between  1  and  1+dl  along  the 
photon's  line  of  flight  is  given  by 

p(l)dl  =  Z^e  dl  (2) 

By  applying  the  basic  principle  to  equation  (2), 

1  -Z  1 

^  *  P(l)  =  /e  Z  ds  =  1-e  (3) 
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solving  for  1 


1  =  j  ln(l-0  (4) 

1  is  then  determined  directly  from  C. 

Equation  (4)  is  often  seen  simplified  to 

1  =  -  I  In  (5) 

^t 

since  l-c  is  distributed  in  the  same  manner  as  C . 

This  example  Illustrates  the  basic  principles,  but  finding  x 
becomes  more  difficult  when  the  Integral  cannot  be  solved  explicitly. 

A  simple  method  to  overcome  an  implicit  integral  is  to  subdivide 
(a,b)  into  Intervals,  storing  accurate  values  of  P(X^)=P^  for  the  end 
points  of  each  subinterval 

X  =  a<x, . . .<x  »  b 

o  n 

Then  using  the  discrete  method  for  determining  the  interval  (X^-1,X^)  on 
which  X  falls.  Interpolation  is  then  used  to  find  the  final  value  of  X. 
The  resulting  equation  is  given  by 


Often  it  is  difficult  or  impossible  to  find  a  closed  form  for  some 
density  functions.  An  alternate  method  called  the  rejection  technique 
may  be  used.  The  rejection  technique  is  often  used  for  sampling  from  a 
density  function  p(x),  a£x£b,  by  employing  two  random  numbers,  ^  and  q. 
Figure  3  illustrates  how  the  technique  uses  the  random  numbers  to  define 
points  distributed  on  the  rectangular  area,  bounded  by  the  lines 
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x=l,  x=b,  y=0,  and  y=l. 


This  technique  consists  of  ’’throwing"  points  (C»n)  uniformly  into 
the  bounded  region  and  rejecting  the  points  lying  above  the  curve  (in 
the  shaded  region).  This  technique  is  best  used  when  the  rejection  area 
is  only  a  small  fraction  of  the  total  enclosed  rectangle. 
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